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Pi  BST8ACT 


This  study  reviews  the  history  ar.d  development  of  marine 
electric  propulsion  drives,  the  types  of  electric  propul¬ 
sion,  and  the  inherent  losses  which  occur  within  the  synch¬ 
ronous  AC  macnir.es  typically  used  for  high-power  propulsion 
systems. 

A  thorough  review  of  the  literature  pertaining  to  heat 
transfer  in  electrical  machinery  is  made.  In  particular, 
the  use  of  liquid  cooling  in  various  flow  configurations, 
including  buoyancy- driven  thermosyphons  and  two-phase  ther- 
aosyphons,  is  analyzed. 

Forced- liqu id  coding  is  feasible,  but  the  required 
rotating  seals  are  a  problem  in  reliability.  Closed-loop 
thermosyphon  cooling  appears  feasible  at  high  rotational 
speeds,  athough  a  secondary  heat  exchange  through  the  shaft 
is  required.  Closed,  two-phase  thermosyphons  and  heat  pipes 
are  also  feasible,  but  require  forced-air  circulation  for 
heat  rejection  to  the  ambient.  since  all  of  these  concepts 
deserve  additional  attention,  areas  for  further  research  and 
development  are  recommended. 
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I.  IHTRODOCTION 


The  aiarine  engineer  and  naval  architect  are  charged  with 
creating  a  design  of  a  vessel  from  a  broad  set  of  require¬ 
ments  for  the  ship  system  prepared  by  the  customer.  This 
requires  a  very  structured  procedure  for  preparing  th® 
design-based  comparisions  and  economic  trade-off  studies  at 
each  and  every  stage  cf  the  development  of  the  design.  The 
end  result  is  a  complete  set  of  design  requirements  from 
which  to  construct  the  final  product  (this  result  should 
meet  the  original  design  requirements!). 

This  structured  procedure  has  the  form  of  the  well-known 
"design  spiral"  [Hef.  1],  which  is  shown  in  Figure  1.1  for 
the  propulsion  machinery.  This  procedure  encompasses  the 
total  ship  design,  including  (but  not  limited  to)  : 

Main  Propulsion  System 
Shaft  horsepower 
Propeller  3PM 

Specific  fuel  consumption  and  bunker  capacity 
Space  and  weight  objectives 
Adaptability  to  ship  configuration 

Auxiliary  Ship  Systems 
Power  and  lighting 

Steam-galley,  deck,  and  heating  systems 
Heating,  ventilating,  and  air  conditioning 


Firefighting  and  ballasting 


Fresh  water 


Hull  Engines  ring  Systems 
Anchor  handling 

Steering  engine  and  bridge  telemetering  control 
Cargo  handling  gear 
Crane  systems 

Electronic  and  Naviqati on  Systems 

Communication,  exterior  and  interior 
Radar 

Loran,  Omega,  etc.,  navigation  aids 

Military  electronics,  sensors,  command  and  control 
systems,  weapons  directors,  tactical  data  systems, 
and  electronic  countermeasures. 

In  order  to  design  the  vessel  within  economic  constraints, 
or  to  use  the  life  cycle  cost  as  the  function  tc  optimize, 
it  is  imperative  to  define  initially  the  devices  used  tc 
provide  the  power  for  the  vessel.  They  must  provide  most  for 
the  least,  in  terms  of  the  fuel  that  they  expend  and  the 
volume  and  weight  that  they  possess. 

The  selection  of  the  main  propulsion  plant  requires 
matching  the  power  of  the  generating  device  with  a  transmis¬ 
sion  (in  most  cases)  ,  the  propulsor,  other  ship  systems  and 

the  hull.  If  one  limits  the  choices  of  propulsors  to  fixed- 
pitch  propellers  and  controllable  and  reversible-pitch 

propellers,  which  are  currently  realistic  for  large  vessels, 
the  number  of  possible  permutations  remains  quite  sizable, 
as  shown  in  Figure  1.2.  The  horsepower  and  RPM  requirements 
cf  modern  warships  require  the  use  of  a  transmission  to 
couple  the  high  RPM  of  an  economical  prime  mover  with  the 
low  RPM  required  at  the  propulsor  for  high  propulsive 
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figure  1.1  Preliminary  Design  Spiral 
for  Propulsion  Machinery. 

efficiency.  This  transmission  may  be  either  mechanical  or 
electrical.  The  mechanical  reduction  gear  has  been 
exploited  in  the  past  to  minimize  volume/space  requirements. 
When  coupled  to  a  reversible  prime  mover,  such  as  a  steam 
turbine,  it  has  '  teen  proven  to  be  quite  acceptable. 
However,  the  electrical  transmission  has  several  attractive 
features,  such  as  ease  of  propulsor  speed  and  direction 
control,  flexibility  in  the  the  number  and  location  of 
controls,  freedom  cf  installation  and  machinery  layout, 
coupling  of  various  prime  movers  independently  to  propulsor 
drives  (as  well  as  ether  auxiliary  applications)  ,  and  th« 
ease  of  auxiliary  thruster  installation  for  maneuvering 
control. 

The  choice  of  the  propulsion  plant  includes  many  impor¬ 
tant  factors,  such  as: 

1.  Reliability 

2.  Maintainability 

3.  Space  and  arrangement  requirements 

4.  Weight  requirements 

5.  Type  of  fuel  required 
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6.  Fuel  consum pticr 

7.  Fractional  power  and  transient  performer.- 
«.  Interrelations  with  -  uxiliaries 

9.  Fever  sing  capability 

10.  Operating  personnel 

11.  Fating  limitations 

12.  Costs 
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Figure  1.2  Propulsion  Plant  Alternatives. 


These  factors  have  gone  through  many  cyclic  changes 
based  on  the  fluctuating  economy.  The  recent  emphasis  on 
fuel  costs  and  availability  has  driven  the  fuel 
consumption/prime  mover  efficiency  to  a  much  higher  level  of 
concern  than  it  was  previously  given.  The  cost  of  acquisi¬ 
tion  cf  new  vessels  has  skyrocketed,  as  well  as  the  costs  of 
maintenance  and  repair.  These  escalating  costs  are  a  chal¬ 
lenge  to  the  engineer  of  today.  The  economic  considerations 
of  the  initial  cost  cf  a  plant,  the  cost  of  maintenance  and 
repair,  and  the  cost  to  operate  a  plant  from  its 


construction  to  its  r etire nest,  also  coined  "lifa  cycle1’ 
cost,  is  the  driving  force  that  controls  the  aesicr.  cf  the 
modern  warship. 

Of  the  list  of  prime  movers  shown  in  Figure  1.2,  the 
life-cycle  cost  f  the  marine  gas  turbine  is  cns  of  the 
lowest.  The  economy  of  operation  is  maximized  within  a 
rather  narrcw  range  cf  RPM's.  This  has  been  tolerated  in 
the  current  employment  of  this  device,  since  the  cost  of 
utilizing  an  electric  transmission  to  accomplish  the  goal  of 
operating  the  prime  mover  within  its  maximum  economy  range 
did  not  offset  the  less  of  efficiency  experienced  by  the  gas 
turbine.  This  was  due  to  the  technology  that  existed  within 
the  field  of  electrical  contrcls  at  the  rime.  The  advent  of 
electronic  frequency  control  of  the  output  of  the  generator, 
independent  of  RPM  cr  poles,  and  the  frequency  control  of 
the  synchronous  motor,  independent  of  its  supply  frequency, 
makes  the  use  of  the  electrical  transmission  much  more 
attractive.  Although  the  initial  cost  will  be  higher,  the 
cost  of  operation,  maintenance,  and  repair  of  the  system 
will  rapidly  offset  the  initial  cost  differential  [Ref.  2], 
[Ref.  3],  and  [Ref.  4].  An  analysis  of  the  economic  and 
engineering  considerations  involved  with  this  mechanical  vs. 
electrical  drive  selection  was  made  by  Eric  Gott  and  S.  0. 
Sver.sson  [Ref.  5]  in  1974;  they  clearly  illustrated  the 
nature  cf  the  problem  and  possible  solutions.  In  *heir 
paper,  it  was  evident  that  the  selection  of  some  of  the 
possibly  more  appropriate  forms  of  electric  drives  was 
discouraged  due  to  the  economics  of  the  choice.  High-power 
frequency-controlling  equipment  was  very  expensive  and 
bulky.  Recent  advances  in  high-power,  solid-state  equipment 
have  now  again  promoted  these  alternatives  into  viable 
contention. 


A.  HISTORIC AL  DEVELOPS  ENT 
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From  the  early  attempts  by  the  Russians  in  the  1890‘s  to 
power  small  passenger  launches,  nr.  oh:  Critic’.  o~  “r • 
experiments  with  small  subm ersibles ,  to  tne  present,  the  use 
of  electric  propulsion  has  mad?  a  slow  evolution.  The  U.S. 
Navy  began  with  the  Jupiter  in  1911.  It  was  a  5500  hp/shaf4- 
vessel  equipped  with  wound  rotor  induction  motors,  which 
fulfilled  a  30-year  lifespan  (terminated  by  surface  warfare 
activity  in  1943,  while  serving  as  an  aircraft  carrier 
(Largely)).  This  led  to  several  major  programs,  summarized 
belo  w : 

*  AC  drive  installed  on  five  battle  ships  cf 
30, 000  +  hp  (1918-1  922) 

*  AC  drive  applied  to  carriers  Saratoga  and  Lexington  of 
200, OOOhp  (1920) 

*  Synchronous  ac  drive  introduced  on  four  Coast  Guard 

*  i 

cutters  (19  19) 

*  Synchronous  ac  drive  in  the  liner  Normandie  at 

160, OOOhp  (  1935) 

*  Synchronous  ac  drive  and  full  dynamic  braking  in 

cutters  (1940) 

In  the  U.S. Navy  Bureau  of  Ships  Manual  (BUSHIPS 
250-660-2)  ,  dated  1  March  1945,  the  text  states  that  the 
Navy*s  inventory  of  ships  that  used  electric  drive  (in  addi¬ 
tion  to  above)  was  ".... about  two  hundred  destroyer  escorts 
and  about  one  hundred  other  ships  including  tankers,  cargo 
transports,  troop  transports,  and  store  ships.”  All  of  these 
applications  were  with  steam  turbines  (which  were  speed- 
adjustable  and  reversible)  and  eventually  this  design  lest 
favor  to  the  less  expensive  mechanical  reduction  gear 
systems  which  were  somewhat  lighter.  The  weight  cf  these 
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origins!  electric  drive  systems  was  due  +  o  oh*  basic 
frequency  relationship  of  REM  to  the  number  cf  poles.  This 
led  to  *  he  basic  configuration  of  two  pole  c^n-retcr? 
(3600RPM)  and  twenty-sight  to  eighty  pole  actors 
(256-9CREM)  .  The  large  number  of  poles  required  in  merer  the 
for  speed  reduction  caused  the  motors  to  be  extremely  large 
ar.d  heavy.  The  advent  of  the  use  of  synchronous  motors 
reduced  the  weight  slightly  since  the  motor  experienced  much 
fewer  slip  losses.  A  good  historical  background  on  the 
early  stages  of  electric  propulsion  is  found  in  the  book 
written  by  Commander  S.M.  Robinson,  OSN  [Ref.  6]. 

Current  applications  in  the  O.S.  Navy  are  shown  in 
Table  I. 


TABLE  I 

O.S.  Navy  Electric  Drive  Ships 


Type 

1  Ships 

SHP 

f  Motors 

AS 

2 

15,000 

1 

AS 

2 

11,520 

8 

TARC 

1 

10,000 

2* 

TAGOS 

1 

1,600 

2** 

*  also  has  4-1200  hp  thrusters 
**  also  has  1-550  hp  thruster 


B.  ELECTRIC  PROPOLSIC1  CONCEPTS 

An  area  of  current  research  in  the  field  of  rotating 
electric  equipment  that  shows  considerable  promise  with 
regard  to  electric  propulsion  is  the  use  of  acyclic  (homo- 
polar)  superconductive  DC  devices.  These  devices  produce 
exceptionally  high  torque  and  horsepower  in  a  very  small 
volume  and  weight.  This  is  of  great  value  in  advanced 
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hul  If  era  propulsion.  A  large  amount  of  research  i.\  r  his 

area  to  determine  configuration,  peak  casualty  jetr  trial , 
maintainability  of  the  required  cryogenic  equipment,  as  well 
as  the  aaaptibility  cf  the  systems  for  this  cryogenic  appa¬ 
ratus  must  be  performed  before  actual  implementation  of  this 
configuration  can  occur  in  a  specific  design.  Thus,  for 
rapid  implementation  in  a  short-term  project  (3-5yrs)  ,  the 
choices  remaining  in  electric  propulsion  to  be  considered 
are  direct-currar.  t  (DC)  and  alternating-current  (AC) 
devices.  The  category  of  AC  devices  is  further  subdivided 
into  inductive  and  synchronous  devices.  All  three  devices 
have  teen  used  for  marine  propulsion  and  each  has  its  cwn 
inherent  advantages  and  disadvantages. 

1 .  CC  Machinery 

The  DC  electric  drive  system  has  been  used  exten¬ 
sively  in  the  past,  primarily  in  low  to  medium  power  ranges 
(1003  -  6000  hp/shaft)  with  extremes  of  400  hp/shaft  and 

IS, 600  hp/shaft  noted.  This  employment  cf  DC  drives 
normally  .utilizes  1-4  prime  movers  powering  one  or  more 
motors  per  shaft  (duplicating  the  system  if  more  than  one 
shaft  exists) . 

Primary  advantages  of  this  system  include: 

1.  Ease  of  Control.  Control  being  effected  by  varying 
the  generator  voltage  through  field  con.trol-a  process 
that  lends  itself  well  to  remote  location(s). 

2.  Boltin le  Control  Stations.  The  ease  of  which  the 
generator  voltage  and  polarity  of  the  motor  connections 
(either  field  or  exciter)  are  varied  allows  multiple 
control  stations  to  be  designed  for  flexible  ship 
control  in  tight  maneuvering  situations. 

3.  Adaptability  to  varing  propeller-hnll  character is- 
£i£S.  The  propeller-hull  characteristics  change 
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dramatically  when  the  vessel  is  towing  or  ice- or  -  ■?.  . 
The  ability  tc  match  the  prime  mover  (rur.nir.  j  at 
cp+,iro"’  c^cor)  ■*•'■  +  '-<o  ore- re  Her  through  the  proDuisior. 
motor  and  still  provide  maximum  torque  is  a  specific 
advantage  of  this  system,  without  excess  engine  or 
electrical  capacity,  through  a  wide  range  of  pr^p  — -e- 
HP  Us. 

Primary  disadvantages  include: 

1.  Limited  Power  Range .  The  top-end  range  of  approxi¬ 
mately  10,000  hp/shaft  limits  the  vessel  size  that  can 
utilize  this  system. 

2.  Generator  S peed  Limitations.  The  maximum  speed 
that  the  generator  can  be  operated,  due  to  problems 
primarily  associated  with  commutation,  is  typically 
limited  to  under  1000  BPS .  This  limits  the  selectior. 
of  the  prime  mover  if  one  does  not  utilize  mechanical 
reduction  gearing  to  match  the  characteristics.  Figure  • 
1.3  shews  the  normal  relationship  of  rpm  vs.  kw. 

3.  Weight  and  Size.  The  weight  of  the  DC  system 
increases  dramatically  with  power  rating  as  shewn  in 
Figure  1.4  (this  figure  shows  weight  in  pounds/SH?)  ; 
note  that  total  weight  becomes  almost  linear  with  SHP 
at  higher  SHP  values.  This  precludes  high  power  appli¬ 
cations  of  conventional  DC  systems  above  10,000 
hp/shaft  (acyclic-homo  polar  devices  are  excluded  from 
this  analysis)  .  The  physical  size  of  the  device  also 
has  this  dramatic  functional  relationship  to  rating  as 
shown  in  Figure  1.5. 

From  the  above  discussion  of  advantages  and  disad¬ 
vantages,  the  consideration  of  DC  electric  propulsion  for 
use  in  a  warship  design  that  would  require  on  the  order  of 
40,000  hp  and  utilize  gas  turbines  in  their  most  efficient 
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Figure  1.5  Variation  of  DC  Generator  Size 

with  Bating. 


inherent  losses  associated  with  these  DC  systems  though  are 
essentially  the  same  as  those  ir.  the  AC  systems.  Figures 
1.1  through  1.5  are  from  "Marine  Engineering",  r.  l. 
Harrington,  editor. 

2 .  AC  machinery 

The  devices  included  in  the  conventional  AC 
machinery  area  fall  into  twc  categories:  induction  and 
synchronous.  The  principal  difference  between  the  twc 
systems  is  in  the  mctor  application.  The  induction  motor 
system  provides  exceptional  torque  at  low  bpm  and  maximum 
slip  (200%).  This  high  torgue  response  at  low  rpm  makes  the 
induction  motor  system  preferable  in  some  applications.  For 
r.aval  aoplicat  ions,  this  system  is  not  as  efficient  as  the 
synchronous  system  by  a  few  critical  percentage  points.  The 


induction  motor  his  losses  that  are  a  function 
incurred  between  the  rotor  and  tha  stator  field. 
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is  r  slated  to  the  RPK  and  load.  The  remaining  losses  ir.  the 
induction  motor  are  the  same  as  the  synchronous  motor  ir. 
general.  Another  critical  difference  between  the  two 
systems  is  the  power  factor  difference.  In  the  induction 
motor,  the  rotor  current  is  induced  from  the  stater  field. 

This  induct!  field  that  is  created  as  a  result  of  this 

curren*  is  the  mechanism  from  which  the  mechanical  rotation 
is  derived.  Since  this  current  is  induced  from  the  stator 
(supplied  current),  it  causes  tha  current  load  seen  by  the 
supplying  device  to  nlag''  the  supplying  voltage  and  there¬ 
fore  cause  the  power  factor  (a  measure  of  the  amount  of  lag 
or  lead)  to  be  less  than  unity.  A  power  factor  of  unity 
indicates  that  the  load  current  and  the  load  voltage  are  in 
phase  and  maximum  real  energy  transfer  is  occurring.  The 
induction  motor  is  normally  operated  with  a  lagging  power 

factor  of  approximately  80?  [Bef.  6]  and  an  overall  effi¬ 

ciency  of  approximately  90%.  In  the  synchronous  system,  the 
power  factor  is  normally  unity  and  is  controllable  by  the 
external  motor  rotor  excitation  which  results  in  the  synch¬ 
ronous  system's  power  factor  being  adjustable  and  a  unity 
power  factor  being  easily  obtained.  Mote  here  also  that  the 
power  factor  in  this  system  can  be  made  to  actually  "lead" 
and  thereby  correct  a  system  "lag”  to  other  loads  that  might 
exist  in  the  fully  integrated  system.  The  advantages  of  AC 
systems  in  general  are: 

1.  High  efficiency.  Losses  in  the  synchronous  system 
are  approximately  6%  (including  excitation  losses)  as 
compared  to  the  induction  motor  system  losses  tha': 


range  from  approximately  8  to  10' 


2.  flexibility  of  Installation.  The  direct  drive 
motor  can  b°  installed  in  the  vessel  to  minimize  shaft 


length  to  the  propeller  (since  physically  smaller  '■hi: 
DC  motors)  and  the  generator  can  be  installed  vi-hcu* 
reduction  gearing  and  operated  directly  coupled  no  th~ 
prime  mover. 

3.  Dual  Ose  of  Pro puls  ion  Power.  The  AC  system  allows 
for  the  power  from  propulsion  units  to  be  used  for 
ether  functions  when  net  required  for  propulsion:  this 
leads  to  the  creation  of  a  fully-integrated  system. 

4.  available  in  large  Power  Ratings.  The  conventional 
AC  systems  can  be  built  in  any  power  rating  required 
with  a  reasonable  upper  limit  of  approximately  60,000 
ho.  The  larger  system's  size  is  being  reduced  by 
raising  the  system  voltage  (which  is  not  limited  by  the 
commutation  requirements  of  DC  systems).  The  system 
voltages  typically  in  use  are  from  2,300  to  about  7,500 
volts.  These  voltages  are  easily  transformed  for  ether 
uses  in  addition  to  propulsion. 

The  selection  of  the  main  propulsion  system  there¬ 
fore  is  not  simple.  The  design  environment  must  be  specified 
and  the  prime  mover  is  selected.  Once  these  parameters  are 
defined,  the  appropriate  selection  of  transmission,  control, 
and  propulsor  can  be  chosen.  For  the  purpose  of  the  propul¬ 
sion  cf  a  vessel  of  approximately  8,000-ton  displacement  and 
the  mission  of  a  normal  warship  of  the  destroyer/cruiser 
variety  using  the  marine  gas  turbine  as  the  prime  mover  and 
life-cycle  costs  the  design  criteria,  the  selection  of  the 
AC  synchronous  drive  (possibly  fully  integrated)  appears  as 
a  very  appropriate  choice.  Some  of  the  design  configura¬ 
tions  possible  within  this  choice,  and  the  associated 
performance  of  each,  are  shown  in  Figure  1.6  as  presented  in 
[Ref.  7].  Additional  possible  configurations  are  shown  in 
Figures  1.7  and  1.8  .  The  end  result  of  the  implementation 
o*  these  technologies  would  result  in  substantial  savings 
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urbine  main  propulsion  in  the  life-cycle  anal; 
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a  comparison  is  shown  in  Figure  1.9,  also 
report  £ Bef .  7 ]. 

Additional  propulsion  plant  arrangements  are  given 
in  [Ref.  5]  and  the  total  number  possible  is  only  limited 
truly  by  the  designer's  imagination  and  creativity  within 
the  confines  of  the  fcullform. 

3.  Hot or/Gene rater  Losses 

The  concept  of  losses  within  the  motor  and  the 
generator  are  comparative  in  magnitude  and  virtually  idan- 
tical  in  nature.  The  losses  will,  therefore,  be  discussed  as 
a  greup  which  is  applicable  to  both  the  motor  and  the  gener¬ 
ator.  The  approach  of  an  energy  balance  is  useful  here  to 
introduce  these  losses: 

| Energy  into  devicel 
|  +  I  on  the  I  = 

I  electrical  side  | 


I  Energy  into  device! 
|  on  the 

I  mechanical  side 


|  Increase  of  I 

|  stored  kinetic/  I 
|  potential  mech.  |  ♦ 
|  energy  in  device  | 


I 

I  Increase  of  I 

|  stored  energy  ! 

|  in  the  elec.&mag.  |  + 
I  fields  of  device  | 


!  Heat  | 
I  gen  erated| 
I  within  | 
I  device  | 


(1.1) 


The  net  quantity  on  the  right-hand  side  of  the  equa¬ 
tion  is  small  and  positive;  the  majority  of  the  energy  into 
the  device  is  transmitted  through  the  device  as  an  output. 
As  an  example,  when  considering  the  motor,  the  energy  into 
the  device  on  the  mechanical  side  is  negative  in  sign  and 
approximately  equal  to  the  electric  energy  into  the  device 
which  is  positive  in  sign.  This  amount  must  be  balanced  by 
the  quantity  on  the  right  side  of  the  equation. 
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Figure  1.7  Integrated  System  Candidate,  Fixed  Frequency 
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Pigare  1.8  Integrated  Systea  Candidate,  Variable  Speed 


Figure  1.9  Comparative  Destroyer  Size  and  Cost. 


The  terms  on  the  right  side  are  all  usually  small  and  may  b* 
zero.  The  first  tern  represents  zhs  energy  associated  wi-.r. 
a  steed  change  ar.d  results  from  the  inertia  of  th«  ■’ 

In  the  steady  state  it  is  zero.  The  second  term  would  result 
from  a  change  in  current  to  the  device  and  the  electric 
field  that  existed  would  be  a  source  of  energy  (the  induc¬ 
tion  effect)  .  The  third  term  is  the  main  point  of  inrerest 
here.  It  is  comprised  of  the  following: 

*  Mechanical  friction  losses 

*  Air  friction  losses  (windage) 

*  Hysteresis  and  eddy-current  losses 

*  Resistance  losses  in  the  armature,  and 

*  Exciter  losses  (in  the  case  of  synchronous  machines). 
The  cbject  for  the  electrical  engineer  is  to  minimize  these 
losses  tc  the  greatest  extent  possible.  Since  it  is  impos¬ 
sible  to  eliminate  all  these  losses,  the  problem  that 
remains  is  the  removal  of  the  resulting  heat  within  the 
device  [Ref.  8].  The  size  and  weight  of  the  device  is 
largely  dependent  on  the  ability  to  remove  this  heat,  since 
the  capacity  of  the  device  is  limited  by  the  insulation  used 
in  the  device.  The  method  of  heat  removal  is  not  limited  to 
a  single  mode,  but  is  a  combination  of  many  modes  of  heat 
transfer  including: 

*  Conduction:  stator  to  casing 

rotor  to  stator  (through  air  gap) 
rctor  to  shaft 
shaft  to  casing 

shaft  to  external  device  (coupling, 
reduction  gear,  etc.) 

*  Convection  (both  forced  and  natural)  : 

casing  to  ambient 

stator  tc  internal  air  (or  gas) 

rotor  to  internal  air  (or  gas) 

shaft  surface  to  internal  air  (or  gas) 
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snarr  surrace  to  external  air 
*  Radiation:  external  surfaces  to  sinks 

internal  surfaces  to  casino 


These  modes  are  r.ot  all  very  controllable.  large 
amounts  cf  research  have  been  devoted  to  minimizing  these 
thermal  resistances.  Efforts  to  improve  convection  heat 
transfer  have  found  that  the  restrictions  on  the  rotor- 
stator  air  gap  have  prevented  exploitation  cf  this  mode  as 
discussed  as  early  as  in  1926  by  G.E.  Luke  [Ref.  9]  and  as 
late  as  in  1979  by  O.N.  Kostikiv,  et  al.  [Ref.  10].  The 
improvement  of  the  conduction  heat  transfer  has  been 
improved  by  minimizing  the  thermal  resistance  cf  the  insu¬ 
lating  compounds  and  the  materials  used  in  the  construction 
cf  the  device,  including  installing  thermally-conduct ive 
materials  in  the  coil  ends. 

The  improvement  of  convection  heat  transfer  is  still 
being  examined  by  using  liquids  and  gases  to  cool  the  rotor 
and  the  stator  of  these  devices.  This  procedure  involves 
piping  a  liquid  through  the  stator  to  remove  heat  by  forced 
convection,  and  utilizing  a  material  with  high  thermal 
conductivity  directly  adjacent  to  the  piping  to  conduct  it 
to  the  piping  where  it  is  removed  by  the  passing  liquid. 
This  method  is  equally  applicable  to  the  rotor,  provided  the 
problem  cf  the  rotating  seals  can  be  overcome  in  order  to 
channel  the  liquid  from  an  external  sink  to  the  moving  rotor 
and  back  out.  The  advantage  of  doing  this  over  standard 
cooling  schemes  is  that  the  capacity  of  a  device  may  be 
increased  due  to  the  lower  internal  temperatures  that  could 
be  maintained  (or  the  physical  size  and  weight  of  the  device 
could  be  reduced  for  an  existing  capacity)  as  shown  in 
Figure  1.10  for  a  40,000  hp,  180  RPM  synchronous  motor.  The 
use  of  a  liquid  in  a  bouyancy  driven  closed-loop,  called  a 
thermosyphon,  is  also  possible.  The  amount  of  heat  transfer 


in  this  method  is  a  a 

*he  liquid  being  used,  the  temperature  rise  allowed,  me 
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the  heat-transfer  areas  then  required  may  be  prohibitively 
large. 

The  two-phase  cooling  method  within  a  closed  loop, 
ir.  the  form  of  a  heat  pipe,  has  also  been  extensively  inves¬ 
tigated.  This  method  has  an  advantage  over  the  conventional, 
ligui d-ccoling  method  in  that  it  eliminates  the  necessity 
for  rotating  seals.  These  seals  can  be  weak  points  in  the 
design  from  a  reliability  stand  point,  especially  if  high 
RPMs  are  involved.  The  heat-pipe  method  typically  involves  a 
confined  fluid  acting  as  a  two-phase  medium  for  heat 
transfer,  exploiting  the  phase  change  as  a  vehicle  for 
substantial  heat  transfer  in  order  to  remove  heat  from  the 
rotor  or  the  stator  cr  both.  This  method  then  transfers  the 
heat  either  to  the  ambient  air  or  to  some  other  gas,  liquid 
cr  solid  ‘  heat  exchanger.  The  configuration  most  suitable 
for  use  within  a  motcr/gen erator  on  the  order  of  40,000  hp 
would  be  either  to  ambient  air  or  to  a  fresh  water  heat 
exchanger  external  to  the  motor  in  the  case  of  stator 
cooling,  or  to  internal  forced  air  through  extended  surfaces 
in  the  case  of  rotor  cooling. 

C.  PROBIEH  DEFIHITIC8 

The  purpose  of  this  thesis  is  to  analyze  various  methods 
of  cooling  motors  and  generators  of  warship  propulsion 
plants  and  to  discuss  the  economic  and  reliability  perspec¬ 
tives  associated  with  each  system.  The  next  two  chapters 
will  discuss  the  theory  and  application  of  these  two 
methods:  liquid  cooling  of  the  rotor,  both  forced  convection 
and  the  thermosyphon,  and  two-phase  cooling  of  the  rotor. 


A  test  apparatus  for  valuation  of  various  ro^a^ir.g 
cooling  schemes  has  teen  designed  and  is  currently  under 
construction  at  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center,  Annapolis,  Maryland.  It  is  sketched  in 
Figure  1.11  and  will  be  used  for  the  analysis  within  this 
study.  The  result s/recom mendations  herein  may  thereby  be 
experimentally  evaluated. 

For  the  purpose  of  this  arbitrary  configuration,  th° 
model  for  analysis  is  as  follows: 

1)  The  device  will  be  a  water-cooled  frame, 
synchronous  AC  generator,  25,118  kVA  a  3,600  RPM. 

2)  A  specific  conduction  bar  shall  be  analyzed  for  a 
typical  load  for  all  cooling  configurations.  It  has  a 
length  of  0.9144m  and  is  0.0116m  sguarewith 

a  4.763mm  diameter  hole  bored  through  the  length 
(herein  referred  to  as  the  tube  length  and  tube 
diameter)  .  The  losses  per  bar  in  this  configuration 
are  approximately  50 W. 
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Figure  1.10  coo?fnC°S§^Imis°n  Various 
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Single-phase  cooling  has  been  extensively  investigated 
with  regard  to  rotating  references.  Work  was  stimulated 
early  by  the  printing  industry,  with  the  need  to  cool 
printing  roll  mills.  The  work  was  further  encouraged  by  the 
gas  turbine  industry  with  the  need  to  cool  high-temperature 
turbine  blades,  which  rotate  at  high  RPM,  and  by  the  elec¬ 
trical  machine  industry,  where  the  need  to  cool  large, 
powerful  generators  became  a  necessity.  Research  in  the 
study  of  flow  and  heat  transfer  in  a  rotating  reference 
system  is  applicable  to  forced-convection  flow,  buoyant  flow 
in  thermcsyphons ,  and  the  single-phase  flows  of  closed  two- 
phase  devices. 

A  theoretical  analysis  of  flow  in  a  heated,  horizontal 
tube  without  rotation  by  Morton  [Ref.  11]  was  completed  in 
1957.  Although  his  analysis  was  limited  to  non-rotating 
flows,  it  gave  considerable  insight  into  buoyar.cy-induced 
secondary  flows.  His  analysis  was  limited  to  laminar 
convection  in  horizontal,  uniformly-heated  tubes  at  low 
Rayleigh  numbers  (based  on  the  temperature  gradient  along 
the  pipe  wall) .  The  work  was  an  extension  of  the  study  by 
Nusselt  who  ignored  the  gravitational  (acceleraticnal) 
effects  and,  therefore,  was  independent  of  orientation. 
This  study  was  a  numerical  solution  to  the  governing  equa¬ 
tions  assuming  small  rates  of  heating  and  a  notable  varia¬ 
tion  of  density  along  the  pipe  length  giving  rise  to 
secondary  flows.  A  similar  study  was  completed  by  Kays 
[Ref.  12]  in  1955.  Morton  studied  this  secondary  flow 

driven  by  the  buoyancy  effects  and  the  resultant  temperature 
distribution  by  successive  approximation  in  a  power  series 
(this  is  similar  to  the  flow  in  a  pipe  rotating  about  a 
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perpendicular  axis).  Seccrd-crder  terms  were 
account  for  the  secondary  flow.  The  basic  equations  and 
d3  v  ~  1  c  r  rn  s  r.  *  i  v.  Appendix  A/  '*?- d  ^ r ?  ^  v  « 

for  many  later  studies.  This  development  is  adapted  tc  the 
rotating  case  by  the  substitution  of  centripetal  accelera¬ 
tion  for  the  gravitational  acceleration  term. 

The  configuration  of  the  cooling  loop  has  three  distinct 
regions  which  must  be  considered:  the  radial  sections  of 

tubing  (oriented  along  the  line  of  the  centripetal  accelera¬ 
tion)  ,  the  entrance  regions  of  the  axial  sections,  and  the 
axial  sections  through  the  conductor  bars  (oriented  parallel 
to  the  axis  of  rotation).  Additionally,  it  must  be  deter¬ 
mined  whether  the  flow  is  laainar  or  turbulent.  The  devel¬ 
opment  will  be  reviewed  chronologically  within  each  topic 
with  sample  calculations  made  for  pertinent  correlations  in 
Appendices  E  and  C  as  noted. 

1.  HISTOHICAL  DEVELOPMENT  OF  ONCE-THROUGH  COOLING 

In  the  196  8  publication,  "Recent  Advances  in  Heat 
Transfer"  (Ref.  13],  Kreith  made  a  survey  of  virtually  all 
the  work  tc  that  date  that  had  transpired  in  the  field  of 
convection  heat  transfer  in  rotating  systems.  His  section 
cn  concentric  cylinders  and  rotating  tubes  summarized  the 
works  mentioned  herein,  although  principally  evaluating  the 
reports  using  gases  as  the  fluid  media.  He  presents  the 
concept  of  utilizing  a  thermosyphon  for  the  coding  of 
rotors,  the  theory  cf  which  is  intimately  related  tc  the 
previous  work  with  single-phase  heat  transfer,  which  will  be 
discussed  later.  Kreith  also  describes  the  three  different 
sections  (radial,  entrance,  axial)  of  the  device  which  this 
study  is  addressing.  These  have  been  analyzed  individually 
both  expermentally  and  theoretically.  The  inlet  region  has 
the  swirl  effect  due  to  the  Coriolis  force,  which  is  a 


function  of  “he  Rossby  number,  or  the  swirl  number  (tr.e 
reciprocal  of  the  Rossby  number)  ;  the  radial  and  axial  tube 
sections  have,  when  a  significant  centrifugal  a.cc^l c  r  a*  ior 
term  is  present,  noticeable  hydrostatic  effects  ar.d,  addi¬ 
tionally,  buoyant  effects  when  a  density  differential 
exists. 

1  •  Padial  Secticr.s 

Extending  his  research  to  a  vertical  tube  rctanina 
about  a  perpendicular  axis,  Morris  [Ref.  14]  confirmed  that 
the  approximate  solutions  must  include  second-order  terms  by 
studying  low  heating  rates  and  using  a  method  similar  to 
Morton,  but  including  tangential  terms  as  well  as  radial 
acceleration  terms.  He  identified  the  three  components  of 
convection: 

1)  Forced  convection  due  to  the  externally-generated 
pressure  gradient, 

2)  Gravitational  buoyancy  in  the  axial  direction,  and 

3)  Rotational  buoyancy  due  to  centrifugal  and  Coriolis 
forces. 

The  relative  magnitude  of  the  terms  due  to  item  3)  is  char¬ 
acterized  by  the  Rossby  number.  This  is  a  ratio  of  th® 
inertia  force  to  the  Coriolis  force.  Although  the  Rossby 
number  exists  in  the  velocity  and  temperature  fields  for 
solutions  through  the  second  order,  it  has  no  effect  on  the 
solutions  for  heat  transfer  or  flow  resistance.  Since  at 
the  nominal  rotation  rates  and  eccentricity  that  exist  in 
rotating  electric  machinery,  the  centrifugal  acceleration  is 
much  much  greater  than  the  gravitational  acceleration,  the 
former  will  dominate. 

In  their  1968  report  on  heat  transfer  in  rotating, 
radial,  circular  pipes  in  the  laminar  region,  Mori  and 
Nakayama  [Ref.  15]  continued  their  research  by  evaluating 
steady,  fully-developed  flow  with  analytical  techniques. 
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Csrr.g  the  tour. i ary- layer  concept,  they  correlated  th *ir 
results  wi*h  experimental  work  by  Trefethen.  They  car* 
correlations  for  the  (Nu/Nu0)  ra-io  as  a  function  of  -he 
angular  velocity,  Reynolds,  Frandtl,  and  rotational  Reynolds 
numbers  for  fluids  with  Prandtl  numbers  above  and  below 
unity  (i.e.,  gases,  liquids,  glycerol,  and  liquid  metals). 
Their  study  confirms  the  stabilizing  effect  of  the  secondary 
flow  cn  the  critical  Reynolds  number  and  notes  that  laminar 
solutions  are  still  valid  when  the  Reynolds  number  is 
considerably  large  (12,600  3  si  =500  rad/s). 

As  a  second  report  on  the  subject  of  rotating  radial 
pipes,  Mori,  et  al.  [Bef.  16]  evaluated  secondary  flow 
effects  due  to  the  Coriolis  force  by  assuming  a  boundary 
layer  alcng  the  wall.  With  experimental  results  and  corre¬ 
lations  from  previous  work,  they  were  able  to  show  that  the 
effect  of  secondary  flow  is  less  in  the  turbulent  region 
than  in  the  laminar  region,  although  still  on  the  order  of  a 
10’S  increase  over  the  value  without  secondary  flow.  The 
correlations  they  presented  agree  with  their  previous  report 
[Bef.  15]  and  with  experimental  results  with  air;  they  are 
presented  as  a  ratio  cf  Nusselt  numbers  based  on  an  exponen¬ 
tial  function  of  the  Reynolds  and  Prandtl  numbers  and  are 
given  for  a  wide  range  of  Prandtl  numbers  including  gases 
and  liquids. 

Ito  and  Nanbu  [Bef.  17]  experimentally  measured 
friction  factors  for  flow  in  rotating,  straight  Fipes  of 
circular  cross-secticn,  and  compared  their  results  to  those 
of  Trefethen  and  Mori  S  Nakayama.  Their  empirical  correla¬ 
tions  for  friction  factor  ratios  for  two  ranges  of  Kt  (where 
Kt  =  J2/Re  and  J=  n  d2/  v  )  are  shown  in  equations  (2.1) 
and  (2.2)  for  the  turbulent  region,  Kt  >  Ktcr.  These  are  as 
follows : 


f/fo- 


0.942  Kto.  282  r  i<  Rt  <500, 


(2.1) 


0.942  Kt 0 • os,  K t >5  00 .  (2.2) 

Eelow  a  Kt  of  1,  the  friction  factor  practically  coincides 
with  tha  stationary  friction  factor  given  by  Blasius: 

f  o  =  °*316  Re-o.as.  <2. 3) 

t 

For  laminar  flow,  which  occurs  whan  Kt  is  below  Ktcr,  the 
following  correlation  applies: 

f/f  =  0.0883  Klo.2S(i  +  11.2  Kl 

U1 

where  Kl  =  J*Re.  Equation 
2.  2*  1 02  <  Kl  <10*  and  (J/Re)<0.5. 
without  rotation  is  given  by: 

f  n  =  64/Re  (2.5) 

U1 

Transition  occurs  at: 

Ktcr  =  1.C7  J1. 23*103  (2.6) 


-0  .  325)  f 


(2.4) 


(2.4)  is  valid  fc: 
The  friction  facto: 


cr 


Re  =1.07  jo-  2  3*1  03  (2.7) 

cr 

for  the  range  28<  j  <2*103.  Calculation  of  the  pressure 
drop  would  be  required  in  the  detailed  calculations  for 
design  of  a  cooling  system,  especially  of  the  clr  •<  !, 
rotat ing-loop  system. 


n 
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Axia^  Sections 

Kuo,  et  al.  [Ref.  18]  analyzed  the  heat-tzar .r f-r 
char acreristics  of  water  flowing  through  a  partially-filled 
pipe  which  was  rotating  about  its  own  axis  and  which  was 
maintained  at  a  constant  heat  flux.  They  also  summarized 
the  previous  pertinent  works  that  were  available  at  the  time 
the  paper  was  written.  The  paper  illustrates  hew  the  tran¬ 
sition  from  gravitational  to  centrifugal  force  is  accompa¬ 
nied  by  instabilities  as  previously  discussed  by  Taylor 
[Ref.  19]  and  [Ref.  20].  The  authors  calculated  a  critical 
Taylor  number  at  which  the  Nusselt  number  changes  its  func¬ 
tional  dependence  on  the  Reynolds  number.  The  paper  also 
shows  how  the  heat  transfer  is  accomplished  within  a  layer 
around  the  periphery  of  the  pipe.  A  correlation  for  this 
film  thickness  is  provided,  and  experimental  results  confirm 
its  validity  and  illustrate  how  the  increasing  depth  of  the 
film  in  the  rotating  pipe,  up  to  completely  full,  does  not 
enhance  the  heat  transfer  beyond  the  critical  depth  value. 
The  range  cf  the  rotational  Reynolds  numbers  used  in  the 
experimental  work  was  between  4,000  and  20,000;  the  enhance¬ 
ment  (Nu/Nuq)  was  approximately  up  to  240*. 

Oliver  [Ref.  21]  analyzed  natural  convection  in  a 
horizontal  tube  with  no  rotation,  following  the  work  of 
Colburn  [Ref.  22],  This  study,  while  pointing  out  several 
problems  with  the  original  research  and  basic  criticisms 
noted  by  Martinelli  [Ref.  23],  combines  the  work  of  many 
authors  to  obtain  a  correlation  (equation  (2.8))  for  the 
Nusselt  number  in  terms  of  the  Graetz,  Grashof,  and  Prandtl 
numbers  and  the  L/D  ratio.  The  paper  also  points  out  that 
the  Grashof  number  must  be  evaluated  using  the  temperature 
gradient  along  the  pipe  axis  instead  of  the  gradient  from 
the  pipe  wall  to  the  fluid  mean  temperature. 

NU  (yw/ub)°‘l*=  1.75{Gzm+5.6*10-*» 
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(Gr  Pr  (L/d)  >  o  .  70}  o.  333. 

m  m 

The  correlation  is  supposed  to  incorporate  the  effsc-s  or 
natural  convection  along  with  forced  convection  by  vector 
addition,  tut  he  nctes  that  the  velocity  profile  of  the 
fluid  flow  is  critical  in  the  analysis  and  than  natural 
convection  in  the  case  of  horizontal  tubes  shc:li  be  ’e- 
pendent  of  L/D  beyond  the  entry  region.  For  this  concept,  a 
correlation  is  also  provided  fcr  the  Nusselr  number  wi*hou- 
the  1/D  term: 

Su  (uw/Mb )  o.i4  =  1.75  {Gzm 

♦  0.0083(GE mPr  )0.7S)o.  333  .  (2.9) 

m  m 

The  study  of  the  Taylor  vortices  was  continued  by 
Pattenden  (  Bef .  24]  using  a  horizontal  (axial)  tube  rotating 
about  its  own  axis.  Due  to  problems  in  the  slip-ring  appa¬ 
ratus  used  in  the  experimental  work,  the  errors  in  the  accu¬ 
rate  measurement  of  the  temperature  precluded  specific 
determination  of  exponents  for  the  Reynolds  number  ir.  the 
correlation  for  the  hea t-transfer -coefficient  equation. 
Qualitative  analysis  of  the  enhancement  due  to  tube  rotation 
was  made,  noting  that  the  axial-flow  effects  were  small 
compared  to  the  rotational  effects.  Taylor  vortices  were 
neither  confirmed,  nor  denied;  the  flow  was  sufficiently 
turbulent  to  mask  the  effects  that  the  vortices  may  have 
had.  Rotation  rat°s  were  varied  to  a  maximum  of  4,000  RPM. 

Mori,  eh  al.  (Ref.  25]  investigated  further  into 
the  effects  of  buoyancy  on  forced-convective  heat  transfer 
in  horizontally-oriented,  heated  tubes  with  laminar  flow. 
His  study  extended  the  work  of  Graetz  and  Nusselt,  and  went 
beyond  the  limitation  of  small  heat  fluxes  such  as  these 
studied  by  Colburn,  McAdams,  Martinelli,  at  al.,  Jackson,  et 
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al.,  a'-n  \’:var.  The  majority  of  the  experimental  work  by 
the  authors  was  c  ?r  formed  wi*h  a  constant  wall  temperature 
and  large  L/0  ratio.  T'.sir  stud/  ir.volvei  a  large  tempera¬ 
ture  difference  between  the  wall  and  the  fluid  such  as  would 
cccur  with  large  heat  flux.  The  change  in  viscosity  with 

temperature  is  approximately  covered  with  the  (y  /y  )Q.  i* 

w  ^ 

term  in  equation  (2.8),  but  the  effect  of  natural  convection 
has  net  been  included  spaci f icaliy.  This  effec*  was  shown 
[Ref.  11]  as  a  function  of  the  PeBa  product;  it  is  appli¬ 
cable  in  the  range  of  PePa  less  than  3,030,  but  within  this 
range  the  effect  of  natural  convection  does  not  exceed 
several  percent.  The  Pavleigh  number  is  a  function  of  the 
fourth  power  of  the  diameter;  therefore,  the  effect  that 
natural  convection  makes  in  small-diameter  channels  is 

small.  This  illustrates  that  in  the  thin  layer  near  the 

wall,  the  temperature  and  velocity  change  abruptly  and  are 
quite  different  from  that  in  the  inner  part  of  the  fluid. 
When  the  RePa  product  is  large,  the  secondary  flow  is  strong 

and  the  velocity  distribution  in  the  axial  direction  is 

utterly  different  from  the  Poiseuille  flow  assumed  in  the 
previous  works.  The  Nusselt  number  is  then  calculated  by: 


Nu  *  2ag/  k  (t  -t  ) 

^  w  m 


(2.  10) 


where  the  stationary  Nusselt  number  for  constant  heat-*flux 


case  is: 


48/11. 


(2. 11) 


The  experiments  were  conducted  with  air  and  fluids  with 
higher  Prandtl  numbers;  they  were  cnly  considering  heating 
in  a  gravitational  field  (not  an  acceleration  field  due  to 
rotation).  In  this  case,  the  layer  near  the  wall  becomes 
thinner  and  the  Nusselt  number  increases  with  the  increasing 


SeSa.  The  local  Nusselt  number  for  the  experimental  corre¬ 
lation  (fer  air)  is  given  by: 


Nu  =  0.61  (SeRa)®*  ^  (1+T.8/ (ReEa)  o*  2) 


(2. 12) 


This  study  was  extended  into  the  turbulent  region,  where  the 
highest  point  of  the  temperature  distribution  was  found  to 
shift  in  the  direction  of  gravity  slightly.  in  this  case, 
the  effect  of  the  secondary  flow  was  not  as  significant  as 
in  the  case  of  laminar  flow.  The  Nusselt  numbers  calculated 
within  this  regime  were  in  good  agreement  with  those  of 
Colburn  [Ref.  22],  where: 


Nu  =  0.0204  Re  0 • 8  . 


(2.  13) 


Mori  and  Nakayama  used  the  theoretical  analysis  of 
their  previous  research  [Ref.  26]  and  [Ref.  27]  on  straight 
pipes  rotating  about  a  parallel  axis  [Ref.  28],  Analyzing 
the  body  forces  driving  the  secondary  flow  caused  by  density 
differences  in  the  centrifugal  field  and  the  Coriolis  force 
with  regard  to  the  flew,  they  used  fundamental  principles  to 
characterize  the  flow  and  temperature  fields.  This  was  done 
assuming  an  effective  secondary  flow  due  to  buoyancy  and 
using  the  rotational  Reynolds  number.  This  would  indicate  a 
rapid  divergence  from  the  results  of  Morris  [Ref.  14]  and 
shows  that  Coriolis  effects  cannot  be  ignored.  The  paper 
gave  correction  factors  as  a  function  of  Reynolds,  rota¬ 
tional  Reynolds,  and  rotational  Rayleigh  numbers,  and 
concluded  that  the  Coriolis  effects  diminish  with  increasing 
eccentricity  (higher  g-fields)  to  a  negligible  value  when 
the  ReRa  product  is  greater  than  10,000,000. 

In  1968,  Nakayama  [  Bef .  29]  further  analyzed  a  hori¬ 
zontal,  straight  pipe  rotating  about  a  parallel  axis.  He 
assumed,  as  previously  [Ref.  15],  that  an  effective 
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secondary  flow  m  fully-d eveloped  conditions  oz  r_ov  : 
temperature  fields  existed,  and  included  body-force  *3ir.s 
explicitly.  The  correlations  he  obtained  for  both  *.h«  ratio 
of  friction  factors  and  the  ratio  of  Nusselt  numbers  are 
valid  over  a  large  range  of  prandtl  numbers  for  both  liquids 
and  gases,  as  well  as  for  a  wide  range  of  Reynolds  and 
Grashcf  numbers.  The  correlation  for  the  Nusselt  number  for 
liquids  in  fully-developed  flow  is: 

Nu  *  Re°*  8  Pr0-*  (0.033  (  Re/  r  2  *  5)  l/30 


(1  +  0.014/(Re/  r2,s)  l/6)) 


(2. 14) 


Inertia  force  is 


where  F  =  Re22/13  (GrrPr 0 •  6) -1  2/1  3 .  The  r  term  is  the 
ratio  of  inertia  force  to  body  force.  Inertia  force  is 
represented  by  Re22/13,  and  the  remaining  terms  in  r 
represent  the  body  force.  The  body  force  is  either  the 

Coriolis  force  in  the  case  of  the  inlet  region,  curved 

regions,  or  the  radial  arms,  or  the  centrifugal  force  in  the 
case  cf  the  axial  section.  The  numerical  result  of  his 

I 

correlations  for  the  model  presented  in  Chapter  I  is 
included  in  Appendix  B  for  both  turbulent  and  laminar 
liquids. 

Siegwarth,  et  al.  (Ref.  30]  further  evaluated  the 
effect  that  bouyancy-induced  secondary  flow  would  have  on 
laminar  flow  in  a  heated,  horizontal  pipe  without  rotation. 
His  work  shewed  clearly  that  secondary  flows  did  exist  due 

to  density  variations  and  that  the  heat  transfer  was 

enhanced  by  them  as  a  function  of  the  Grashof  and  Prandtl 
numbers.  This  effect  is  somewhat  negated  in  the  presence  of 
extreme  acceleration  as  that  in  the  periphery  of  a  rotor, 
but  remains  nonetheless,  provided  the  turbulence  is  below 
the  critical  Reynolds  value. 


Kith  no  assumptions  as  to  the  structure  cf  the  i 2c v 
and  temperature  fields  to  simplify  the  governing  equations, 

of  the  rotating  cylindrical  tube  (Ref.  31]  with  laminar, 
fully-developed  flow.  Data  were  compared  to  the  theoretical 
results  for  the  case  of  air,  water,  and  glycerol.  The 
correlation  obtained  is  in  terms  of  the 

Rayleigh-Reyr.olds-Pr andtl  product  as  is  given  by: 

Nu  =  Nu0  (0.262  (RaRePr)  0,1  73}  .  (2.15) 

The  value  for  the  test  model  is  calculated  and  presented  ir. 
Appendix  B. 

Stephenson  again  studied  this  parallel,  rotating- 
pipe,  heat-transfer  problem  for  fully-developed  turbulent- 
flow.  He  compared  his  experimental  results  and  correlation 
with  the  earlier  work  of  Morris  and  Woods  [Ref.  32]  and  they 
compared  favorably  within  the  turbulent  entrance  region.  He 
noted  that  the  rotational  buoyancy  was  not  a  strong  factor 
in  the  secondary- flow  in  the  fully-developed  region. 
Because  of  this,  his  results,  when  compared  against 
Nakayama's  results  assuming  a  strong  effect  [Ref.  29],  did 
not  compare  well.  His  correlation  is  listed  below: 

Nu  =  0.0071  RoO.68  JO. 023  (2.  16) 

and  is  compared  with  previous  correlations  for  the  test 
model  in  Appendix  B. 

3 •  Entrance  Regions 

In  1969,  the  Institution  of  Mechanical  Engineers, 
Thermodynamics  and  Fluid  Mechanics  Group,  sponsored  a  sympo¬ 
sium,  on  the  subject  of  heat  transfer  and  fluid  flow  in 
electrical  machines.  Included  in  the  presentations  was  a 


Faper  by  David  sen  [  Bef .  33],  which  described  a  fall- sic-., 
gene  rater-rotor  test  rig  fer  evaluating  the  heat-transfer 
characteristics  of  hydrogen  cooling.  His  testing  showed 
that  enhancement  in  the  turbulent  flow  regimes  was  net  as 
high  as  theoretically  predicted.  It  was  hypothesized  that 
♦he  reason  was  the  inability  of  the  hydraulic-diameter 
concept  to  accurately  accour^  for  the  flow  conditions  in  the 
non-circular  ducts.  The  tests  were  able  to  give  a  good 
turbulent  correlation  with  theoretical  prediction  when  the 
cooling  scheme  was  modified  by  shortening  the  axial  path. 
This  effect,  although  not  explained,  could  have  been  due  to 
the  Coriolis  swirl  effect  in  the  inlet  region  enhancing  the 
heat  transfer  above  that  to  be  obtained  in  the  fully- 
develcped  region. 

The  effects  of  the  entry  length,  especially  with  low 
Prandtl  numbers  such  as  those  of  gases,  are  very  difficult 
to  eliminate  and  the  Coriolis  effects  on  the  heat  transfer 
may  be  noted  even  at  high  Reynolds  numbers.  As  the 
Eeynolds-Rayleigh  product  increases  above  about  1,000,000, 
there  is  a  tendency  for  the  amount  of  enhancement  to 
diminish.  This  is  attributed  to  the  turbulent  effects  over- 
iding  the  secondary  flow  enhancement.  The  effect  of  rota¬ 
tion  on  heat  transfer  is  to  enhance  it  in  the  entrance 
regions  by  the  secondary  flow  due  to  the  Coriolis  effects 
and  in  the  fully-developed,  horizontal-tube  regions  by  buoy¬ 
ancy  effects.  A  similar  survey  of  the  technology  was  ntad<=> 
by  Petukhov  and  Polyakov  in  the  USSR  with  much  the  same 
results  [Ref.  34],  The  heat-transfer  problem  in  the 
entrance  region  of  a  tube,  where  the  Coriolis  terms  domi¬ 
nate,  is  addressed  by  Horris  and  Woods  [Ref.  32]. 
Correlations  are  presented  for  air  and  are  applicable  for 
ether  gases  as  well.  These  correlations,  for  both  the 
laminar  and  the  turbulent  cases,  are  presented  as  functions 
of  the  product  of  the  Reynolds  number  and  the 


49 


-V*— . 


T 


.■» T'T -7  '• . r-'.  *"■  ST  •*.'  .*  1 


>1.  ,  '  II  ,.  .Jl_  Ul_  ^  |.W 


rr-.  »  1  .  u-,.  v 


i  ■  .  *  J  >  ■  .  ■  | 


r!~  t  a  t  dor.  al  —  5  sy  n  o  id  s  number.  The  Morris  and  Woods'  rarer 
notes  that  further  work  is  required  in  the  case  of  liquids, 
although  a  similar  approach  is  valid,  Mori  arc  ;.r-.k  ayama 
also  presented  a  survey  of  the  stat a-of-the-art  technology 
regarding  the  heat-transfer  characteristics  of  rotating 
pipes  and  ducts  [Ref.  35].  The  questions  of  the  entrance 
effect  and  its  length  were  emphasized  as  being  generally 
unanswered,  but  they  gave  a  general  guideline  for  deter¬ 
mining  the  extent  of  the  entrance  region  as  approximately  20 
times  the  diameter  of  the  pipe.  They  also  verified  correla¬ 
tions  cf  previous  works  for  helium  and  water. 

4.  Combinations  cf  Ra dial.  Horizontal.  arjd  Entrance 

Sections 

Also  presented  in  the  1969  Institution  of  Mechanical 
Engineers  Symposium  was  a  study  by  Lambrecht  [Ref.  36], 
which  discussed  the  problems  associated  with  water  cooling 
of  rotors.  His  paper  summarized  the  theoretical  considera¬ 
tions,  as  well  as  reviewed  his  previous  work  and  the  work  of 
fellow  German  researchers  (Neidhoeffer  and  Ingenieurs)  .  He 
noted  the  superior  cooling  properties  of  water  and  the  size 
and  weight  reduction,  along  with  the  improved  efficiency  of 
water-cooled  machines.  He  listed  values  for  the  optimum 
cooling  duct  size  based  on  the  heat-tr  ansf  er  and 
electrical-loss  characteristics  for  various  gases.  The 
method  used  for  optimum  duct  calculation  was  based  on  the 
principle  that  the  pressure  drop  in  the  duct  decreases  with 
increasing  duct  size  and  electrical  losses  increase  with 
decreasing  conductor  area  (increasing  duct  size)  .  A  suit¬ 
able  calculation  of  this  type  is  necessary  in  any  final 
design  for  cooling  of  these  devices  by  any  method  devised 
for  either  the  single-  or  the  two-phase  approach. 

In  1970,  Sakamoto  and  Pukui  measured  heat-  and 
mass-transfer  coefficients  for  air  and  oil,  specifically  for 


us?  iti  c co ling  rotating  «l°ctric  machines  [Her.  37].  mini 
that  the  most  effective  and  economic  method  to  increase  nil 
power  output  of  electrical  machines  was  tc  improve  th'~ 
cooling  of  *he  insulating  marerial  (especially  in  th«  rotor 
conductors  and  drums) ,  a  geometry  similar  to  a  typical  rotor 
was  fabricated  and  suitable  coolants  were  tested.  The 
convective  heat  transfer  was  characterized  as  a  function  of 
the  L/d,  Re ,  3rr,  Pr  variables,  with  the  H/d  parameter  held 
constant.  When  the  shape  factor  of  eccentricity  H/d  is 
sufficiently  large,  the  centrifugal  acceleration  may  be 
assumed  to  be  a  uniform  acceleration  (  2H)  acting  on  the 
tuba  axis,  which  is  similar  to  the  force  gravity  exerts  on  a 
horizontal  tube.  with  a  lew  Reynolds  number,  the  empirical 
correlation  was  in  the  form: 


Nu  =  Nu8  (1+C  (Rar)  m/Gz)  n. 


(2.  17) 


where  the  values  of  the  constants  were: 


C  =  0.03  m  =  0.75  n  =  1/3 

and  where  Nu  =  1 . 86  (RePrd/L)  33  for  large  L/d  values, 
and  Nu  =  2.67  Gz0*33  for  small  L/d  values. 

This  was  valid  within  +10*  to  -17%  for  Reynolds  numbers 
between  162  and  2,700  and  rotational  Rayleigh  numbers  to 
20,000,000.  It  was  also  found  that  the  length  of  the 
entrance  region  seemed  to  have  no  affect  on  the  correlation, 
which  wculd  imply  that  a  secondary  flow  existed  in  both 
configurations.  An  interesting  point  in  this  paper  is  that 
the  rotational  Grashcf  number  (hence  rhe  rotational  Rayleigh 
number)  was  calculated  using  the  temperature  difference 
between  the  wall  and  the  fluid  instead  of  the  axial  tempera¬ 
ture  gradient. 

Woods  and  Morris  investigated  laminar  flow  in  the 
rotor  windings  of  directly -cooled,  electrical  machines  in 
197u  [Ref.  38]  and  noted  that  if  water  was  used  as  the 


m  =  0.75 


coolant  m  lieu  of  hydrogen  (which  is  currently  widely 
used),  an  efficiency  improvement  of  >0.5  'h  could  be  octaiiVid. 
Given  that  the  0.5*  cx  a  large,  say  1,000  MW,  generator  was 
the  only  consideration  made  (and  neglecting  the  size  and 
weight  reduction  and  the  extended  operating  life  possible 
due  to  the  cooler  operation) ,  they  proceeded  to  analyze  the 
fundamentals  of  the  problem.  They  surveyed  the  previous 
works  directly  applicable  to  this  area  of  coding; 


1)  Morris  [Bef.  14]  used  a  series-expansion  technique 
(valid  fcr  low  rotation  rate  and  low  heating  rates)  ,  and 

2)  Mori  and  Nakayama  [Bef.  28]  assumed  a  secondary  flow 
(which  was  claimed  to  be  valid  fcr  high  rotation  rates)  ,  and 
used  an  integral-type  analysis. 

Woods  and  Morris  claimed  both  were  inadequate  due  to  the 
restrictions  imposed  ty  the  nature  of  their  solutions.  They 
attempted  a  versatile  and  "exact"  solution  for  laminar  flow 
in  the  fully-developed  region  by  solving  the  governing  equa¬ 
tions  wi+h  a  numerical  procedure.  Numerical  solutions  were 
presented  for  both  the  friction  factors  and  'the  Nusselt 
numbers.  Their  analytical  values  were  compared  with  experi¬ 
mental  results  and  to  correlations  of  previous  works. 
Eiscr epancies  were  explained  and  the  difficulty  in  obtaining 
fully-developed  conditions  in  the  experiments  of  their  work 
and  previous  works  was  emphasized.  The  axial  density- 
variation  affects  were  noted,  also,  as  a  potential  for  error 
in  their  analysis  at  high  values  of  rotational  F.ayleigh 
numbers. 

Nakayama  and  Fuzicka  evaluated  the  generator 
problem,  where  reasonably  large  radii  of  rotors  (approxi¬ 
mately  1  m)  and  high  BPM  (3600)  create  high  centrifugal 
acceleration  effects  [Ref.  39].  This  acceleration  causes 
the  centrifugal  buoyancy  term  to  be  significant.  It  also 


emphasizes  the  Coriolis  term  ir.  the  radial  rube  s*c3?r.t?  and 
curved  inlet  regions.  They  presented  correiat ions  for  rhs 
radio  of  friction  factors  and  the  radio  of  Nusscld  mb»;sf 
arriving  at  fried! cr.-f acoo r  raoics  for  each  of  the  three 
regions  (radial,  curved  inlet,  and  axial),  and  compared  the 
results  to  experimental  data.  Their  correlations  are  listed 


fcslo v : 

f/f0  =  2.2  Ro-1.33  (radial)  (2.13) 

f/fo  =  1.5  Ro~o*3  (coil  end)  (2.19) 

f/fo  =  14.0  Ro-o.8  (axial)  (2.20) 

Nu  =  Re° •  3  Pr0**  {0.014  (Re/Ro2*5)  0.12*)  (radial)  (2.21) 


Their  report  confirms  the  results  of  Nakayama's  correlation 
(equation  (2.14))  for  the  axial  sections. 

The  application  of  the  technology  of  rotating 
cooling  schemes  is  net  limited  to  the  printing  industry  and 
electric  machine  area;  extensive  research  has  been  devoted 
to  this  field  of  study  by  the  gas  turbine  industry.  The 
advisory  Group  for  Aerospace  Research  and  Development 
(AGARE)  held  a  meeting  in  September,  1977,  at  Ankara, 
Turkey,  to  discuss  High  Temperature  Problems  in  Gas  Turbine 
Engines.  Numerous  papers  were  presented  on  the  subject 

including  a  paper  by  W.D.  Morris  on  flow  and  heat  transfer 
in  rotating  coolant  channels  [Bef.  40].  He  used  a  selection 
of  experimental  results  to  illustrate  the  influence  of 
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rotation  on  heat  transfer,  and  demonstrated  that  Coriolis 
and  centripetal  inertial  effects  significantly  altar  tin 
heat -transfer  characteristics  compared  to  the  stationary 
case.  Using  the  governing  equations,  he  illustrated  hew  the 
Navier-Stokes  equations  are  modified  to  include  the  acceler¬ 
ation  effects  of  the  Coriolis  forces  and  the  centripetal 
forces.  For  constant  properties,  these  equations  then  show 
how  at  distances  sufficiently  downstream  from  the  entry 
region,  where  the  axial  gradients  of  the  velocity  are  negli¬ 
gible,  the  elimination  of  the  pressure-gradient  terms  from 
the  radial  and  tangential  momentum  equations  causes  the 
Coriolis  terms  to  vanish  identically  as  a  source  for  the 
creation  of  secondary  cross-stream  flow;  see  also  [Ref.  14]. 

In  the  entry  region,  the  Coriolis  terms  interact 
with  the  developing  axial  velocity  and  create  secondary 
flows  perpendicular  to  the  axial  direction,  even  with  the 
constant  property  condition.  Finally,  in  the  fully- 

developed  region,  the  affect  of  the  rotation  manifests 
itself  only  as  a  cross-stream  pressure  distribution  if  buoy¬ 
ancy  terms  are  neglected.  If  the  density  is  allowed  to  vary 
with  temperature,  the  centripetal  acceleration  terms  need 
only  be  included  as  the  Coriolis  terms  do  not  affect  this 
action.  This  cross-stream,  buoyancy-induced  flow  gives  rise 
to  greatly-enhanced  beat  transfer  and  an  attendant  increase 
in  flow  resistance.  The  effects  of  the  eccentricity,  with 
regard  tc  the  Coriolis  terms,  has  little  effect,  provided 
that  the  ratio  of  the  radius  (H)  to  the  diameter  of  the  tube 
(d)  is  greater  than  5  (only  very  small  radius  of  rotation 
will  have  any  noticeable  effects). 

Marto  [Ref.  41]  lists  the  papers  presented  at  the 
14th  Symposium  of  the  International  Centre  for  Heat  and  Mass 
Transfer  (ICHflT)  held  in  Dubrovnik,  Yugoslavia,  in  1982. 
Papers  of  particular  interest  to  this  study  included  a  paper 
by  Jchnscn  and  Morris  (Ref.  42],  wherein  the  concept  of 
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csn'-.ripetal-acceleration-induced  secondary  flow  was  orcvsr. 
to  be  cf  little  importance,  and  its  influence  was  shewn  to 
be  tctallv  hydrostatic.  Their  work  also  confirm  ad  the 
effect  cf  Coriolis  acceleration  with  regard  to  secondary 
flows,  especially  in  sections  with  relatively  small  L/d 
ratios.  No  quantitative  recommendations  were  presented. 

5 •  Flo w  Transit  ions 

Mori  [Ref.  25]  states,  when  considering  turbulent 
flow,  that  it  is  not  necessary  to  consider  the  influence  of 
buoyancy  on  heat  transfer.  The  secondary  flow  was  found  to 
suppress  the  turbulence  level  when  the  turbulence  at  the 
inlet  region  was  high  and  an  empirical  formula  for  rhe  crit¬ 
ical  Reynolds  number,  in  terms  of  the  Reynolds  number  and 
Rayleigh  number  product,  was  calculated: 

Re  =  12  8  (Re  Ra)  o-  (2.  22) 
cr 

• 

Contrary  to  this,  when  the  turbulence  at  the  inlet  region 
was  low,  the  critical  Reynolds  number  was  higher  (7,700  vs. 
2,000),  and  heating  decreased  the  critcal  value.  The  net 
result  was  that  when  the  ReRa  value  was  high,  the  secondary 
flow  caused  by  buoyancy  makes  the  critical  value  of  Reynolds 
number  tend  toward  the  same  value,  whether  the  turbulence  in 
the  inlet  is  high  or  low.  in  the  case  of  low  turbulence  in 
the  inlet,  the  value  for  the  critical  Reynolds  number  was 
given  by: 

Re  =  Re  /( 1+0.  14  ReRa*10-5)  .  (2.23) 
cr  cr9  '  '  ' 

The  graph  shown  as  Figure  6 2  in  the  Krieth  article  [Ref.  13] 
is  shewn  below  as  Figure  2.1.  This  graph  shows  the  heat- 
transfer  as  a  function  of  both  the  RPM  and  the  flow-through 
Reynolds  number.  It  is  easy  to  see  the  flow's  transitional 


dependence  on  the  Reynolds  number  as  a  function  of  ppm.  ror 
the  experimental  model  being  evaluated,  where  P.  e=t ,  uo  C  a.t  - 
Rar=7,4'?o(t  the  transition  is  delated  to  P.e=11,000  for  -h~ 
:orrelation  that  assumes  high  turbulence  at  the  entrance. 
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Figure  2.1.  Effect  of  Rotational  Speed 
in  a  Rotating  Pipe. 


The  transition  to  turbulent  flow  is  suppressed  by 
rotation  and  has  been  qualitatively  noted  by  numerous 
authors.  An  exact  correlation  for  transition  to  turbulent 
flow  in  a  rotating- reference  frame  with  heating  is  r.ct  known 
to  exist.  From  the  available  literature,  it  is  clear  that 
the  transition  follows  a  path  from  the  laminar  correlations 


to  the  turbulent  correlations 


This  path. 


the  rotating  reference,  has  a  ’'dip”  a 
in  a  plot  of  the  Musselt  number  against 
Figure  2.2,  whereas  in  the  non-rota: 
occurs. 


in  the  case  of. 
shewn  qua lit at iv ti¬ 
the  Reynolds  number, 
ng  case,  r.c  "dip" 


Re^_  >  Re  2  >  Re^  >  Re  4 

Figure  2.2  Transitional  Characteristics  of 
Botatmg  Systems.  , 


6 .  Summary 

Eased  upon  the  information  reviewed  in  the  litera¬ 
ture,  it  may  be  concluded  that  the  three  regions  within  the 
motor/generator  application  (shown  in  Figure  1.11)  have  very 
different  heat-transfer  characteristics. 

1)  In  the  radial  section  cf  the  coolant  path,  the  Coriolis 
forces  tend  to  dominate  the  heat  transfer  and  friction. 
This  effect  is  to  enhance  the  heat  transfer  and  increase  the 
pressure  drop  in  these  sections.  The  heat-transfer  enhance¬ 
ment  is  substantial  and  could  be  effectively  exploited  for 
the  heat  rejection  from  the  device  by  axial  fcrced-air 
convection  of  the  substantially-hollow,  synchronous  rotor. 


2)  Ir.  the  coil-end  regions  (both  the  radial  cent  -  ~k-  -m 
points  ar.d  t ha  coil-end  loops)  the  Coriolis  fore*-  is  stmo 
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result  is  enhanced  heat-transfer  and  fricticr. -factor 
effects,  but  to  a  lesser  degree  than  in  the  radial  sections. 

3)  In  the  winding  bars,  the  centripetal  forces  dominate; 
these  forces  result  in  hydrostatic  effects  which  limit  flows 
that  would  be  affected  by  the  Coriolis  forces  ar.d  result  in 
some  degree  of  secondary  flow  due  to  buoyancy.  The  amount 
cf  secondary  flow  is  questionable  and  is  somewhat  affected 
by  the  degree  of  turbulence  and  the  amount  enhancement  is 
lass  specific. 

The  enhancement  in  the  laminar-flow  heat  transfer 
has  been  reported  by  numerous  authors.  This  can  be  seen  in 
Figure  2.3.  Both  the  Nakayama  [Ref.  29]  and  the 
Woods-Morris  [Ref.  31]  correlations  for  Nusselt  number  fall 
between  the  classical  laminar  value  of  48/11  and  the  turbu¬ 
lent  value  cf  the  Dittus-Boelter  correlation.  They  have  the 
same  trend  Throughout  the  range  of  RPM  shown  and  differ  by 
only  a  few  percentage  points.  Apparently,  the  effect  of 
rotation  in  laminar  flow  enhances  the  heat  transfer  because 
cf  intense  secondary  flows  due  to  centripetal  forces.  An 
interesting  point  is  that,  at  the  higher  RPM  values,  the 
correlations  of  Nakayama,  Woods-Morris,  Stephenson,  and 
Dittus-Bcetler  all  are  exceptionally  close  together. 

The  heat-transfer  predictions  for  turbulent,  fully- 
developed,  parallel  flows  are  shown  in  Figure  2.4.  This 
figure  illustrates  that  the  resultant  heat  transfer  is  some¬ 
where  in  the  vicinity  of  the  classical  turbulent  correlation 
cf  Dittus-Bcelter  for  the  current  model. 


tor  (turbulent) 


Figure  2.3  Variation  of  Laminar  Nusselt 
Correlations  with  RPH. 


The  hea*  transfer  is  either  somewhat  enhanced  (on  *  he  cr de  ~ 
cf  30?)  in  the  lower  RPM  ranges  (Nakayama  [Ref.  29  ]  )  ,  or 
somewhat  suppressed  ty  the  high  centripetal  accsl  =  raticr. 
forces  as  shown  by  the  correlation  of  Stephenson  [Ref.  43]. 

The  enhancement  of  the  heat  transfer  in  the  Nakayama 
correlation  [Ref.  29]  in  the  lower  RPM  range  is  assumed  to 
be  due  to  the  residual  Coriolis  affects  from  the  entrance 
regions  remaining  longer  in  the  lower  centripetal  accelera¬ 
tion  field. 

B.  HISTORICAL  DE7EL0PHEHT  OF  THEBHOS YPHON  COOLING 

The  thermosyphon  is  a  device  that  creates  the  motion  of 
the  working  fluid  solely  by  buoyancy  forces.  It  has  a 
region  where  heat  is  supplied  to  the  working  fluid  and  a 
region  where  heat  is  rejected.  The  less  buoyant,  cooler 
portion  cf  the  fluid  moves  in  the  direction  of  the  accelera¬ 
tion  vector,  while  the  opposite  is  occurring  for  the  warmer 
portion.  The  thermcsyphon  also  has  a  third  region,  the 
coupling  region,  with  problems  of  fluid  shear,  mixing  and 
entrainment.  Figure  2.5  illustrates  the  conceptual  arrange¬ 
ments  of  the  three  basic  thermosyphon  systems. 

Kreith  [Ref.  13]  mentioned  in  his  review  that  the  ther- 
mosyphon  (i.e.,  a  liquid -filled  container)  had  a  great 
potential  for  cooling  of  electric  rotating  machinery.  In 
motors  and  generators,  the  use  of  a  thermosyphon  would 
preclude  the  necessity  of  using  rotating  seals,  which  have 
marginal  reliability  in  applications  of  high  RPM  and  volu¬ 
metric  flow.  In  1973,  Japikse  [Ref.  44]  gave  a  thorough 
review  of  the  literature  regarding  thermosyphon  technology. 
Beginning  with  a  reference  to  the  Davies  and  Morris  paper 
[Ref.  45],  which  will  be  discussed  in  detail  later,  Japikse 
categorizes  thermosyphons  as  follows: 


1)  the  nature  of  boundaries  (either  opened  or  do?**  to 
mass  flow)  , 

2)  the  regime  of  heat  transfer  (purely  natural  convec¬ 
tion  cr  mixed  natural  and  forced  convection), 

3)  the  number  and  types  of  phases  present,  and 

4)  the  nature  of  the  body  force  present  (gravitational 
or  rotational)  . 

These  broad  categories  are  then  restriced  to 
prescribed  circulating  fluid  system  driven  by  thermal  buoy¬ 
ancy  forces.” 

The  use  of  the  closed  thermosyphon  in  rotating  refer¬ 
ences,  such  as  in  gas  turbine  cooling,  has  been  studied  in 
detail.  In  1965,  Bayley  and  Lock  (Kef.  46]  experimentally 
investigated  the  performance  of  closed  thermosyphons.  They 
utilized  the  theoretical  development  of  the  open  thermosy¬ 
phon,  modified  the  theory,  and  experimentally  verified  their 
results.  It  was  shewn  that  the  critical  operating  parame¬ 
ters  include:  , 

1)  Length-to-Diameter  Ratio.  This  ratio  controls  the 
characteristics  of  the  flow  and  determines  the  extent  of  the 
coupling  region. 

2)  Heated-Length-to-Cooled-Length  Ratio.  As  this  ratio 
tends  toward  zero,  the  analysis  is  identical  to  the  open 
thermesyphon  correlations.  This  restricts  the  use  of  the 
analysis  for  the  closed  thermosyphon  to  Lh/Lc  ratio  to  be 
not  much  greater  than  one. 

3)  Coupling  Region.  This  is  the  region  that  compli¬ 
cates  the  analysis  and  causes  the  Prandtl  number  to  affect 
the  heat  transfer.  It  exists  in  three  distinct  modes; 
conduction,  convection,  and  mixing.  These  are  regarded  as 
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ideals  and  actual  coupling  usually  consists  of  a  combination 
of  all  thre°. 

Thermosyphon  technology  was  also  applied  to  transfer:?: 
cooling,  nuclear-reactor  cooling,  heat-exchanger  tins,  heme 
and  industrial  furnaces,  cryogenic  cool-down  apparatus, 
steam  tubes  for  bakers'  ovens,  internal  combustion  engine 
cooling,  and  environmental  control  in  the  space  vehicle 
program,  as  well  as  many  others  mentioned  in  the  Japikse 
article. 

In  1971,  Bayley  and  Martin  [Bef.  47]  also  reviewed  the 
state-cf-the-art  technology,  with  particular  emphasis  on 
gas-turbine  applications.  They  studied  both  the  open  and 
closed  thermosyphon  systems.  The  use  of  an  open  thermosy¬ 
phon  system  in  a  rotating  reference,  such  as  in  gas  turbine 
cooling,  as  well  as  in  the  cooling  of  electric,  rotating- 
machinery  area,  has  problems  regarding  the  fluid  selection 
and  its  saturation  temperature  and  pressure.  In  the  open, 
rotating  reference  system,  the  region  that  receives  the  heat 
is  under  high  acceleration  forces  with  high  resultant  pres¬ 
sures.  As  the  natural  flow  due  to  buoyancy  occurs,  the 
heated  fluid  moves  from  a  high  pressure  to  a  relatively  low 
pressure  area  where  a  phase  change  to  vapor  may  occur, 
possibly  even  explosively,  blocking  the  flow.  They  intro¬ 
duced  the  idea  of  using  two-phase  cooling  within  the  closed 
ther mosvphor.  as  a  way  of  exploiting  this  phenomenon,  by 
using  the  high  heat  fluxes  associated  with  the  latent  heats 
of  evaporation  and  condensation  and  the  lower  temperature 
gradients  associated  with  the  phase  change  processes.  This 
also  has  an  added  benefit  of  reduced  weight  over  the 
single-phase  liquid  system.  They  noted  that  further 
research  in  the  area  of  the  Coriolis  force  effects  on  the 
heat-transfer  chact eristics  in  a  rotating  reference  needed 
to  be  accomplished . 


Applying  the  technology  of  the  ther mosyphon  to  electric 
machines,  Morris  and  Davies  [Hef.  45]  adapted  the  principles 
of  the  closed  thermcsyphon  to  a  closed  rotating  loop  tha* 
could  be  placed  into  the  rotor  assembly  of  an  electrical 
device.  This  configuration  avoids  the  entry  choking  and 
nixing  that  occurs  in  the  typical  closed  thermcsyphon  and 
generally  avoids  the  complex  mid-tube  exchange  process,  as 
well  as  the  adverse  core-boundary  layer  interaction.  The 
strong  acceleration  force  present  in  the  periphery  of  the 
rotor  of  electric  rotating  machinery,  coupled  with  the 
heating  due  to  the  electrical  losses,  induces  stong  loop 
circulation.  This  heat  can  then  be  rejected  either  to  the 
shaft  of  the  rotor  acting  as  a  heat  exchanger  itself,  tc  the 
air  passing  through  the  hollow  rotor  region  and  the  air 
passing  through  the  end-bell  region  of  the  rotor  ends,  cr  to 
both.  The  shaft,  acting  as  a  heat-exchanger  itself,  can  be 
cooled  with  an  internal,  rotating  heat-pipe  (or  two-phase 
thermosyphon)  which  could  transfer  .the  heat  axially  to 
another  heat  sink  outside  of  the  motor/generator  casing.  As 
an  alternative,  the  shaft  could  be  internally  coded  by 
forced  convection  of  water  (cr  some  other  fluid)  to  an 
external  heat  exchanger  through  the  ncn-coupling  end  of  the 
shaft  by  the  use  of  an  axial  rotating  seal.  The  reliability 
and  the  casualty-control  aspects  of  the  axial  rotating  seal 
make  it  preferable  tc  the  radial  rotating  seal  mentioned 
earlier. 

Theoretical  analysis  of  the  loop-flow  has  used  simple, 
one-dimensicnal  force  and  energy  balances,  and  need  net  be 
repeated  here.  The  closed  loop  is  shown  in  Figure  2.6  as 
conceptually  applied  to  a  rotor.  Of  course,  numerous  loops 
would  be  designed  into  the  rotor  to  carry  away  the  heat 
generated.  Incorporation  of  a  one-way  valve  was  suggested 
as  being  necessary  tc  insure  one-directional  flow  created  by 
the  mean  temperature  difference  and  the  pressure  differences 


Figure  2.6  Rotating  Loop  Thermosyphon  for  a  Rotor 


in  the  lints  of  the  circuit.  Experimental  testing  was 
conducted  in  a  test  rig  similar  to  the  rig  proposed  for  the 
evaluation  of  the  systems  studied  herein  and  now  under 
construction  at  DTNSBDC.  A  modest  range  of  rotation  rates 
was  employed  in  the  Hcrris  and  Davies  test  rig,  typical  of 
most  meter  applications.  The  details  of  their  test  section 
were  very  similar  tc  the  model  used  for  analysis  herein: 
25.4  mm  (on*  inch)  diameter  copper  rod  through  which  a  6.35 
mm  (0.25  inch)  hole  was  bored.  The  experimental  results 
were  placed  in  the  form: 


i'u  -  -  (Gi.ljP.Oi  P<.  f  Ac)  • 


i  -■  i 

,  i .  -  -«  j 


The  ratio  of  angular  acceleration  to  gra  virarional  accelera¬ 
tion  (Ac)  was  the  authors’  way  of  incluling  the  rotational 
effects  (previous  authors  have  used  the  rotational  Grashof 
number).  The  final  correlation  was: 


Nu  =  0.1505  aco*73s  Re2.»s  pr/Gr 


(2.25) 


for  the  range  of  103<  NuGr/Pr  <106. 

The  critical  consideration  in  the  analysis  of  the 
closed-loop  thermosyphon  is  the  temperature  rise  in  the 
heated  limb  (conductor  bar)  .  The  model  described  in  Chapter 
I  was  analyzed  by  using  the  governing  relationships  from  the 
Davies  and  Morris  [Ref.  45]  paper  iterated  against  the  Woods 
and  Morris  [Ref.  31  ]  correlation  for  laminar,  fully- 
developed  flow  in  a  heated  circular  tube  rotating  about  a 
parallel  axis.  It  can  be  seen  from  Figure  2.7  that  the  heat 
transfer  increases  with  RPM  (the  higher  acceleration 

TABLE  II 

Temperature  Rise  in  the  Clos9d-Loop  Thermosyphon 


1000 

3000 

3600 

5000 


Tem£.  Rise  (  C) * 
209.  0 

93.5 
41.0 

35.5 
27.  0 


*  for  a  cold-side  temperature  of  45  C. 

effect)  .  For  the  conditions  on  the  model  which  constructed 
Figure  2.7,  the  values  of  temperature  rise  is  shown  in  Table 
II  for  various  RPM's.  This  clearly  shows  that  this  method 


motor  device  (i.e. ,  lew  RPS)  ;  on  the  other  hand,  it  coal! 
be  very  reasonable  for  the  generator.  The  friction 
was  assumed  to  be  64/Re  and  no  correction  was  mads  for  th= 
rotation  in  the  friction- factor  calculation;  the  minor 
losses  were  neglected.  If  the  friction  factor  is  increased 
due  tc  the  rotation,  as  in  equation  (2.4),  and  minor  losses 
(laminar)  are  included,  the  temperature  rise  will  be 
increased  slightly.  The  trend  will  be  the  same  and  the 
temperature  rise  at  the  higher  RPM's  will  still  be  accep¬ 
table.  The  actual  design  of  the  closed-loop  thermosyphon 
will  require  an  optimization  of  the  heat -transfer  calcula¬ 
tions  and  the  friction- fact  or  calculations.  This  optimiza¬ 
tion  may  indicate  that  a  larger,  single,  clossd-thermcsyphor. 
loop  for  a  conductor  bar  group  is  preferable  over  individual 
conductor  loops, 

C.  ADVANTAGES/DIS AD VANTAGES  OF  LIQUID  COOLING 

The  main  objectives  of  cooling  are  to  remove  the  heat 
created  by  the  electrical  losses  and  improve  the  operating 
efficiency,  extend  the  lifetime  of  the  insulation,  and 
reduce  the  overall  weight  and  size  of  the  device.  The 

feature  cf  liquid  coding  (vs.  gas  cooling)  is  the  superior 
capacity  of  liquids,  especially  water,  to  remove  heat  with  a 
small  rise  in  temperature.  .  The  principle  disadvantage  of 
using  liquids,  or  gaseous  fluids  other  than  air,  is  the 
design  of  the  piping,  coupling,  and  secondary  heat  exchanger 
required  to  support  the  device  cooled  by  the  fluid.  Since 
one  of  the  main  objectives  is  to  reduce  weight  and  size,  the 
fluid  chosen  must  be  able  to  remove  enough  heat  per  unit 
mass  to  accomplish  this  goal  within  the  device. 
Additionally,  since  the  application  under  evaluation  is 
marine  combatants,  the  fluid  should  be  non-toxic,  ar.d  th® 


Figure  2.7  Theraosyphon  Analysis  for  Nusselt  number 


system  should  be  durable  enough  to  withstand  damage  an? 
maintain  operational  safely.  The  popular  fluid  tc  date  ir. 
land-based  ir.st allaf ions  has  been  hydrogen  which  is  cl*arlv 
an  unacceptable  choice  for  the  marine  combatant  due  tc  its 
explosive  nature  and  since  the  size  reduction  obtained  is 
not  sufficient  to  warrant  system  development.  Of  the 
remaining  fluids,  the  obvious  choice  for  its  high  heat 
capacity,  availability,  and  safety,  is  water.  The  use  of 
water  has  been  tasted  both  in  stators  and,  on  a  limited 
scale,  ir.  rotors. 

The  principle  problem  to  overcome  in  the  use  of  water  is 
in  the  rotary  seals  needed  to  direct  the  fluid  intc/cut  of 
tha  machine.  Additionally,  when  the  system  inside  the 
device  is  analyzed,  the  ability  of  a  design  to  withstand 
damage  is  questionable.  The  design  would  be  critically 
dependent  on  the  fluid  within  the  rotor  for  operation,  even 
in  a  casualty  mode.  The  supply  of  cooling  water  to  a  stator 
would  no*  be  as  critical,  since  any  casualty  could  be 
repaired  easily  with  shipboard  damage-control  equipment,  and 
emergency  cooling  of  a  stator  could  be  accomplished  by 
external  means  due  to  the  thermal  conductivity  of  the  stator 
housing.  For  the  rotor,  this  leaves  a  choice;  either  not 
utilize  the  enhanced  cooling  and  suffer  the  weight  and  size 
penalties,  or  develop  a  clo  ^d  system  for  cooling  the  rotor. 
The  closed  system  for  the  rotor  is  not  necessarily  isolated. 
A  closed-loop,  thermcsyphon  could  be  placed  in  the  rotor, 
transferring  the  heat  from  the  rotor  to  the  rotor  shaft. 
This  loop  could  have  additional  heat  transfer  to  air  passing 
through  the  core  of  the  rotor.  Further,  the  heat  trans¬ 
ferred  to  the  shaft  could  he  removed  to  the  ambient  air  or 
an  external  heat  exchanger  utilizing  the  non-coupling  end  of 
the  motor  by  the  use  cf  circulating  water  through  the  shaft, 
or  by  a  rotating  heat  pipe  as  seen  in  [Bef.  48]. 


An  additional  method  for  cooling  the  rotor,  which  ha? 
been  studied  by  numerous  authors,  is  the  use  cf  hear  cipes 
or  closed,  two-phase  thermo  syphons.  In  this  case,  the  hea* 
is  transferred  from  the  rotor  to  air  circulated  through  the 
end-bell  areas  of  the  device.  This  method  cf  rotcr  cooling 
greatly  enhances  the  survivability  of  the  device  in  a 
damage-ccntrol  sense.  Each  heat  oipe  is  independent  of  the 
others  and  damage  to  a  few,  unless  all  in  the  same  group  of 
adjacent  conductor  bars,  would  not  japardize  the  device; 
indeed,  even  if  all  were  in  the  same  group,  limited  opera¬ 
tion  could  continue. 


III.  TWO-PB  ASE  FLUID  COOLING 


The  concept  of  utilizing  the  latent  heat  of  vaporization 
to  transfer  heat  is  widely  used  in  heat  exchangers  every¬ 
where#  from  household  and  industrial  boilers  and  space 
heaters  to  marine  propulsion  units.  Air-conditioning  equip¬ 
ment  and  heat-pump  units  also  employ  this  idea  to  advanta¬ 
geously  utilize  the  low  temperature  of  vaporization  of 
fluorc-carbcn  compounds  to  remove  the  heat  from  the  ambient 
and  release  this  through  condensation  under  pressure,  thus 
transferring  the  heat  in  the  desired  direction  at  the  cost 
of  the  energy  to  pump  and  pressurize  the  gas  to  the  point  of 
condensation.  In  each  of  these  two-phase  applications,  the 
transport  of  the  working  fluid  may  require  a  pump,  or  it  may 
be  accomplished  by  the  buoyant  force  that  exists  in  an 
acceleration  field,  such  as  that  due  to  gravity.  This  is 
the  principle  that  is  employed  in  the  closed  two-phase  ther¬ 
mosyphon.  This  device  is  similar,  to  the  classical  "heat 
pipe"  and  the  only  difference  is  that  the  "pumping"  method 
in  the  heat  pipe  uses  capillary  action  instead  of  accelera¬ 
tion  forces  to  pump  the  fluid.  Thus,  a  "heat  pipe"  used  in 
a  rotating  reference  utilizing  the  acceleration  forces  to 
transport  the  liquid  from  the  condenser  section  to  the 
evaporator  section  (in  lieu  of  a  wick)  is,  technically,  a 
two-phase  ther mosyphcn.  Figure  3.1  illustrates  the  opera¬ 
tion  of  the  closed  two-phase  thermosyphon.  The  evaporator 
end  must  be  in  the  direction  of  the  acceleration  field  for 
the  liquid  to  be  transported  from  the  condenser.  Figure  3.2 
illustrates  the  operation  of  the  heat  pipe,  which  utilizes 
the  capillary  action  of  the  wicking  material  to  transport 
the  liguid  from  the  condenser  to  the  evaporator.  The 
remainder  of  the  chapter  shall  use  these  definitions  for  a 
discussion  cf  both  devices. 
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V-V 


Acceleration 


Field 


Figure  3.1  Typical  Two-Phase  Thermosyphon. 


Acceleration  Field  Not  Required 


Figure  3.2  Typical  Heat  Pipe. 


1.  HISTORICAL  DE7ELCPBEIT 

The  development  cf  the  heat  pipe  (which  will  apply 
equally  to  the  development  required  for  the  two-phase  ther- 
mosyphon)  was  begun  by  Gaugler  [Ref.  49]  in  1944  and  by 
Trefethen  (Ref.  50]  in  1962.  Trefethen  was  working  on 

cooling  methods  for  spacecraft  (zero-g  environment)  for 
General  Electric  Company  and  discovered  that  capillary 
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pumping  can  be  a  very  valuable  area  for  further  development. 
Working  independently,  Grover,  at  ai.  [Ref.  51]  r-inv-r. ned 
this  concept;  he  widely  published  the  first  applications  cf 
the  device  and  gave  it  its  name.  Kraus,  et  al.  (Ref.  52] 
have  reported  this  historical  development  and  have  given 
sample  calculations  for  the  design  of  a  heat  pipe.  Chi 
[Ref.  53]  also  reviews  the  theoretical  development  of  heat 
pipes,  and  provides  the  design  procedures,  including  sample 
calculations,  for  their  utilization. 

B.  CONVENTIONAL ,  CLOSED  TWO-PHASE  THEHMOSTPHON  ANALYSIS 

As  previously  mentioned,  the  appealing  factors  associ¬ 
ated  with  the  two-phase  utilization  include;  high  heat 
fluxes  associated  with  phase  changes,  lower  temperature 
gradients  associated  with  these  processes,  and  reduced 
weight  of  the  two-phase  system  over  the  single-phase  liquid 
system. 

Research  by  Cohen  and  Bayley  [Bef.  54]  (referred  to  and 
discussed  by  Japikse  [Bef.  44]  )  found  that  the  amount  of 
liquid  filling  the  system  functionally  affected  to  the  heat 
transfer  in  both  rotating  and  static  tests.  They  found  that 
heat  transfer  increased  as  the  percent  liquid  in  the  system 
was  increased  to  approximately  1.5%  by  volume,  then 
decreased  to  an  intermediate  value,  and  finally  increased  to 
approximately  the  same  value  as  in  the  1.5%  case.  This  was 
related  to  the  following  process:  in  the  low  filling  situ¬ 
ation,  condensate  returns  to  the  evaporator  section  and 
forms  a  thin  film  on  the  walls.  It  is  in  this  film  that  the 
heat  transfer  occurs  and  the  phase  change  takes  place.  In 
the  case  of  the  completely-filled  evaporator,  a  liquid  pool 
region  develops,  and  a  condition  of  nucleate  boiling  exists. 
If  the  filling  is  insufficient,  with  regard  to  the  heat  flux 
being  transported,  the  pool/film  will  not  continue  to  the 
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end  of  -the  heated  section  and  "dry -out”  will  occur.  This 
•'dry-cut”  will  assaae  a  drop-wise  instead  of  film-wise  phase 
change.  The  geometry  of  the  regions  in  which  the  films  are 
likely  ic  form,  influences  the  results. 

Lee  and  Mital  (Ref.  55]  conducted  experiments  with  an 
electrically-heated,  water-cooled  thermosyphon  using  water 
and  freon  as  the  test  fluids  and  varied  the  filling  quan¬ 
tity,  L,  /L  .  pressure  (or  T  .  ) ,  and  heat  flux  (L,  /L  is  the 
ratio  cf  the  length  of  the  heated  section  to  the  cooled 
section)  .  The  result  of  the  filling  quantity  on  heat 
transfer  was  the  same  as  that  reported  by  Cohen  and  Bayley: 
increasing  heat  transfer  with  filling  to  a  point  and  then 
decreasing  beyond  that  value  to  an  intermediate  value  and 
increasing  to  the  case  of  the  completely-filled  evaporator. 
The  effect  of  decreasing  was  to  increase  the  heat 

transfer  within  the  range  0.8  to  2.0  and  the  advantage  of 
larger  condenser  area  was  evident.  The  heat-t ransf er  coef¬ 
ficient  was  found  to  increase  significantly  with  increasing 
mean  pressure  due  to  at  least  three  factors: 


since  the  mass  flew  for  a  given  heat  flux  is  nearly 
constant  and  density  of  the  liquid  increases  with 
pressure,  a  lower  pressure  drop  (and  lower  AT)  is 
necessary  fer  the  same  heat  flux; 


for  a  larger  ?sat #  ?sat  var^6S  more  rapidly  with 

T  for  the  fluids  considered,  hence  requiring  smaller 


AT's  at  higher  pressures;  and 


for  lower  pressure  drops,  a  more  favorable  force  balance 
exists  on  the  condensate  film  permitting  a  faster  liquid 


return. 


Water  was  found  to  give  heat  fluxes  superior  to  those  of 
freon,  for  for  the  same  at,  due  to  the  larger  values  of 
latent  heat  of  vaporation,  and  thermal  conductivity. 


V.v.V 


Lee  and  Mital  [Ref.  55]  also  considered  the  ar.aly-ical 
problem  of  predicting  the  maximum  heat-transfer  rate  for  a 
laminar  film,  constant-wall-temperature  condenser  and  a 
constant-heat- f 1 ux  evaporator.  Neglecting  the  forces  due  to 
vaper  pressure  drop  and  momentum  changes,  and  using  a  force 
balance  on  the  falling  film,  balancing  the  effects  of 
gravity  and  fluid  shear,  they  related  the  mass  flow  rate  to 
the  heat  flux.  Osing  a  local  energy  balance  and  an  overall 
energy  balance  between  the  condensing  section  and  the  evapo¬ 
rator  section,  the  following  relations  for  the  heat  trans¬ 
ferred  (g)  and  the  saturation  temperature  (T  )  were 

sat 

obtained : 

q  *  {p2R3h  g/2u  Lhgc}c  <3.u 

fg  n  c 

and 

q  =  {k(T  -T  )/(RL  /L  )}C/D  '  (3.2) 

sc  n  c 

where  C  =  1/8  -  1/2  y2  *  3/8  y4  -  1/2  y4  In  y  , 

D  *  y4ln  y  (1/2  In  y2  -1/2)  ♦  1/8  y4  +y2(1/2  In  y  -1/4)  =1/8 

and  y  =  1  -  d/R.  The  simplifying  assumptions  (forces  due 
to  vaper  shear  and  momentum  changes  were  neglected  in  the 
force  balance)  cause  an  error  (of  as  much  as  a  factor  cf  2) 
to  occur,  but  qualitative  behavior  is  correct  as  far  as 
Lh/Lc ,  Tsatand  working  fluid  are  concerned.  This  develop¬ 
ment  is  also  included  in  the  survey  by  Japikse  [Ref.  44], 
This  heat  transfer  is  graphically  presented  in  the  Lee  and 
Mital  paper  in  their  Figure  12,  shown  here  as  Figure  3.3. 
This  clearly  shows  that  the  quantity  of  heat  that  can  be 
removed  by  a  two  phase  water  device  is  large,  even  at  a  low 
saturation  temperature.  The  report  also  indicates  that  for 
water  the  increasing  ratio  of  heated  length  to  cooled 
length,  L+,  decreases  the  maximum  heat  transfer.  The  heat 
transfer  is  also  very  sensitive  to  the  operating  pressure; 
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Figure  3.3  Comparison  of  the  Experimental  Results  with  the 
Analytical  Prediction-from  Lee  5  Mital  paper. 


the  heat -transfer  coefficient  increases  appreciably  with  the 
mean  operating  pressure.  The  temperature  drop  increases 
with  the  pressure  drop  along  the  length  of  the  thermosyphon 
tube  (and  the  saturation  pressure  gradient).  As  the  pres¬ 
sure  increases,  the  vapor  specific  volume  decreases, 
resulting  in  a  decrease  ir.  the  pressure  drop;  the  mass  flow 
rate  of  the  vapor  is  essentially  constant  for  transferring  a 
given  heat  input  rate.  In  the  previous  papers,  the 
vertically-oriented  thermosyphon  was  considered  and  the 
current  irodel  (for  seme  of  the  orientations)  requires  the 
device  to  be  oriented  in  the  horizontal  direction  (perpendi¬ 
cular  to  the  applied  acceleration  field)  .  Although  the  Lee 
and  Mital  paper  was  considering  a  vertically-oriented,  two- 
phase  thermosyphon,  the  general  results  are  the  same  for  the 
horizontal  tube. 
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C.  LISITS  OP  OPER&TICfl 

When  considering  the  design  of  a  two-phase  thermcsv then, 
or  a  heat  pipe,  limitations  to  the  heat  transfer  must  be 
considered;  four  of  these  limits  are  common  to  both  devices 
(sonic,  entrainment,  toiling,  and  condensing).  A  qualita¬ 
tive  comparison  of  these  heat-transfer  limitations  as  a 
function  of  the  saturation  temperature  is  given  in  [Sef.  48] 
and  is  shown  in  Figure  3.4. 


Figure  3.4  Operating  Limits  of  Rotating  Heat  Pipes. 

The  sonic  limit  and  the  boiling  limit  have  been  readily 
analysed.  The  entrainment  and  the  condensation  limits  are 
not  as  well-known  and  little,  if  any,  literature  describing 
these  phenomena  exist.  This  is  an  area  in  which  further 
research  is  needed. 

1)  Sonic  limit;  the  vapor  flow  in  a  two-phase  closed 
limited  to  the  sonic  velocity  at  the  operating 


system  is 


The  sonic 


pressure  (also  known  as  the  "choking"  limit). 
liai4-  is  represented  by: 


Q 

max 


A  V  , 

v  s  fg 


(3.3) 


where  Vg  is  the  sonic  velocity  of  the  vapor.  This  limiting 

velocity  may  be  determined  experimentally,  computed  by  the 

relationship  V  =((3p/3p.)  }  l/2 ,  or  approximated  bv  tha 
s  s 

perfect-gas  relationship  {YHT  }  »/z. 

sat. 

2)  Entrainment  limit :  the  interfacial  shear  between  the 
liquid  and  the  vapor  will  held  the  liquid  (which  is  flowing 
in  the  opposite  direction)  back  and  starve  the  evaporator 
section  cf  liquid.  This  counter-current  flow,  when  the 
relative  velocity  is  large,  causes  the  interface  tc  become 
unstable  which  results  in  waves  at  the  interface.  As  the 
vapor  velocity  continues  tc  grow,  droplets  of  liquid  are 
formed  at  the  liquid  surface  as  the  shear  force  exceeds  the 
surface-tension  force.  The  formation  of  these  droplets  and 
their  subsequent  entrainment  in  the  vapor  stream  causes  the 
partial  or  total  stoppage  of  the  flow  (dry-out) .  This 
phenomenon  is  generally  governed  by  the  Weber  number  (the 
ratio  cf  the  inertial  force  to  the  liquid  surface-tension 
fcrce) .  The  Froude  number  is  also  used  to  characterize  the 
phenomenon  cf  the  drep-wise  entrainment  of  the  liquid  in  the 
vapor  stream.  In  the  application  considered  herein,  the 
formation  cf  the  waves  is  considered  unlikely  due  to  the 
extreme  acceleration  field  present.  Since  the  flow  is 
counter-current  and  is  in  the  presence  of  a  high  accelera¬ 
tion  field  the  entrainment  limitation  thus  will  reduce  to 
the  "hold-up"  limit.  It  will  remain  stratified  until 
hold-up  occurs,  resulting  in  evaporator  dry-out. 
Experimental  determination  of  the  exact  correlation  for  the 
hold-up,  or  entrainment,  limit  is  required  for  a  horizontal 
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tuba  rotating  around  a  parallel  axis,  thereby  creating  a 
perpendicular  acceleration  field  many  times  greater  than 
gravity.  The  entrainment  limit  is  generally  given  by: 

^max=^fgpvAvVE»  (3.«) 

where  is  when  the  Weber  number  equals  unity,  the  Froud? 
number  equals  unity,  or  is  experimentally  found.  In  another 
approach,  Jaster  and  Kosky  [Bef.  56]  used  the  ratio  of  the 
axial  shear  force  of  the  vapor  flow  on  the  liquid  surface  to 
the  gravitational  (accelert ional)  body  force  upon  the  liquid 
(defined  as  "F")  in  order  to  arrive  at  an  appropriate 
"hold-up"  criterion.  They  correlated  experimental  data  as  a 
function  of  this  "F"  value,  as  the  flow  they  were  studying 
transitioned  from  stratified  to  annular.  The  result  of 
their  experiments  was  that  the  flow  was  stratified  for  "F" 
values  below  five  (5)  and  was  f ully-annular  above  twenty- 
nine  (29).  The  value  of  "F"  set  equal  to  five  (5)  was 
established  as  the  criterion  for  the  transition  from  strati¬ 
fied  to'  annular  for  the  co-current  flow  case.  The  same  type 
of  analogy  was  developed  by  Collier  and  Wallis  [Ref.  57]  who 
balanced  the  inertia  force  and  the  acceleration  force  to 
scale  stratification  effects  where  their  criterion  was  given 
as  j*  (volumetric  flux)  equal  to  0.25.  Both  of  these  refer¬ 
ences  point  out  the  need  for  further  research  in  this  area. 

3)  Boiling  li*it:  Again,  as  discussed  in  the  first 
section  of  this  chapter  in  the  general  discussion  on  two- 
phase  thermosyphons,  the  boiling  limitation  in  the  current 
application  is  due  tc  the  creation  of-  profuse  nucleaticn  at 
the  evaporator  section  and  a  resultant  vapor  film  on  the 
evaporator  surface  which  insulates  the  evaporator  wall, 
resulting  in  dry-out  and  overheating.  The  boiling  limit  is 
given  by  the  Zu  ber-Kutateladze  prediction  as  discussed  in 
(Bef.  48]  and  is  given  as: 


Q  =0 .  1 3 A  o  1  / 2 h  {  ,Q2R(p  -o  )  a  }  »/ 

max  3'  v  fg  1  v 


(3.5) 


»v.«re  A  is  the  area  in  t ha  evaDorator  section  for  liasiio 
B 

heatir.q. 


4)  Condenser  liyit:  The  capacity  of  the  condenser  is 
dependent  cr  many  parameters,  such  as  the  geometry,  the 
weeding  fluid,  the  orientation  of  the  vessel  with  respect  to 
the  acceleration  field,  and  the  operating  conditions.  The 
orientation  of  the  vessel,  for  the  model  considered,  is  that 
of  a  horizontal  tube  (perpendicular  to  the  acceleration 
field)  with  internal  condensation.  The  condensing  limit  is 
the  capacity  to  condense  the  vapor  on  the  inside  surface  of 
the  condenser  section,  and  this  is  shown  in  Figure  3.4.  No 
analysis  on  the  outer  surface  is  presented  in  this  thesis. 
Instead,  the  condenser  limit  is  based  on  the  capacity  of  the 
inside  wall  exposed  to  condensation,  using  an  arbitrary,  but 
reasonable,  wall  temperature.  (Generally,  the  largest 
thermal  resistance  cccurs  between  the  ambient  and  the 
condenser  cuter  surface.  Therefore,  heat  transfer  is 
limited  by  this  resistance  unless  the  outside  surface  area 
is  substantially  increased  by  the  use  of  fins.)  The 
condensing  limit  is  given  by  an  equation  from  [Ref.  58]  by 
Collier  for  a  horizontal  tube  with  internal  condensation  in 
a  gravitational  field  as: 


Q  =A  ATF  {(p  (P  -p  )  gh  k  *)/ 
max  v  c  l  1  v  fg  f 


(3.6) 


For  the  rotating  reference  being  considered,  the  "g"  term  in 
equation  (3.6)  is  replaced  by  the  centripetal  acceleration 
of  q  2R .  The  factor  "F  "  in  this  equation,  allows  for  the 
fact  that  the  rate  of  condensation  on  the  statified  layer  of 


liquid  is  r.egligable.  The  value  of  "F  "  depend?  or.  rh^ 
angle  p  and  is  tabulated  in  [Bef.  58].  Its  value  ranges 
from  zero,  when  the  tube  is  full  (no  surface  remains  for 
condensation  to  occur),  no  0.725  ,  when  the  rube  ampry.  The 
angle  $  is  shown  in  Figure  3.5. 


Figure  3.5  Laminar  Condensation  within  a  Horizontal 

Tube  (from  collier)  . 


1.  Computation  of  Heat-Transfer  Limits 

The  calculation  of  these  limits  is  an  important  step 
in  the  design  process. 

The  calculations  required  to  analyze  the  configura¬ 
tion  of  the  current  model  are  not  all  well  known.  The 
following  appear  to  be  the  most  applicable  to  the  model 
being  evaluated: 

- The  sonic  limit  calculation  is  from  basic  principles  as 

shown  by  equation  (3.3)  by  allowing  the  vapor  velocity  to 
increase  to  the  Mach  value  of  unity. 

- The  boiling  limit  is  given  by  the  Zuber- Kutateladze 

prediction  as  discussed  in  [Bef.  48]  and  is  given  by  equa¬ 
tion  (3.5) 


- The  condensing  limit  is  given  by  an  equation  from 

(Ref.  58]  ty  Collier  in  equation  (3.6). 

The  literature  reveals  no  correlation  yet  for  the  entrain¬ 
ment  limit  in  the  heated,  horizontal  tube  rotating  about  a 
parallel  axis  with  counter- current  flow.  Since  the  majority 
of  analyses  previously  conducted  for  this  type  of  flow  have 
shown  that  the  entrainment  limit  is  typically  below  the 
sonic  and  above  the  condensing  limit,  the  correlations 
previously  mentioned  have  been  plotted  against  the  amount  of 
fluid  that  is  charged  in  the  tube  (percentage  fill)  for  the 
three  reasonably  known  values  (sonic,  boiling,  and 
condensing  limits)  in  Figure  3.7.  Also  shown  is  the 
entrainment  limit  computed  from  the  Jaster  and  Kosky 
[Ref.  56]  correlation  developed  for  co-current  flow. 

As  was  previously  noted,  the  condensing  limit  is  thQ 
factor  that  is  the  most  limiting  in  the  design  of  the  ther- 
mcsyphon.  The  values  of  Qmax  are  increased  to  a  maximum 
based  on  the  length  of  the  condenser,  the  percentage  fill 
and  the  A  t  available  between  the  constant  temperature 
(assumed)  wall  and  the  saturation  temperature  of  the  vapor 
condensing  on  it.  As  already  mentioned,  the  resistance  to 
heat-transfer  afforded  by  the  condenser  outer  walls  to  the 
ambient  may  be  quite  significant.  The  ability  to  design  an 
adequate  fin  arrangement  and  supply  cooling  air  in  suffi¬ 
cient  quantities,  while  minimizing  the  windage  losses, 
itself  is  a  substantial  problem.  From  the  data  obtained  in 
this  rough  model,  there  are  some  self-evident  points  to  be 
considered : 

- The  length  of  the  condenser,  from  a  heat-transfer  point 

of  view,  determines  the  heat-removal  potential, 

- The  amount  of  liquid  fill  can  be  optimized. 
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- The  temperature  at  which  the  walls  of  the  cor.d-.nser  ar ~ 

maintained  and  the  temperature  at  which  the  then mosyp her.  is 
operated  are  critical  as  this  Temperature  difference  is  -h 
driving  force  in  the  heat  transfer. 

for  the  model  being  evaluated,  and  neglecting  (as 
the  limiting  value)  the  entrainment  limit,  the  closed, 
•►wo-phase  t her mosv phen,  operated  at  a  saturation  temperature 
at,  or  above,  100°C  and  with  the  condenser  section  held  at  a 
wall  temperature  of  50°  C  with  rotation  at  a  0.381  m  radius 
and  3,600  RFM  can  transfer  the  required  50  W  of  power  with 
the  condensing  limit  being  the  limiting  value.  Figure  3.7 
is  based  on  a  heat  load  of  50  W  to  be  equally  divided 

between  condenser  sections  in  both  end-bell  regions  (i.e. , 
approximately  25  W  each  per  condenser)  .  For  the  3600  RPM 
model  the  percentage  fill  may  increase  to  approximately  40*, 
and  still  maintain  the  required  heat  flux.  The  effect  that 
the  RFM  has  with  regard  to  the  temperature  rise  in  the 
heated  section,  makes  the  choice  of  the  closed-loop  thermo¬ 
syphon  pcor  for  the  motor  application,  but  quite  suitable 

for  the  generator.  Further  background  on  the  theory  of 
two-phase  flow  is  available  in  publications  by  Wallis 

[Hef.  59]  and  by  Collier  5  Wallis  in  [Ref.  57].  The  feasi¬ 
bility  of  the  design  has  been  shown  by  Corman  and  McLaughlin 
[Ref.  60]  and  by  Groll,  et  al.  [Ref.  61]  through  experi¬ 
ments  with  heat-pipe  cooling  of  electric  motors  in  the 
scheme  primarily  recommended,  although  the  difficulty  of  the 
liquid  transport  at  lew  rotation  rates  was  noted.  A  diagram 
of  the  cooling  scheme  such  as  recommended  and  tested  by 
Corman  8  McLaughlin  and  Groll,  et  al. ,  is  shown  in  Figure 
3.6  (reproduced  from  [Bef.  60]  )  . 


84 


Entraining  Limit 


Figure  3.  7  Heat-Transfer  Limits  in  the  Closed 
Two-Phase  Tnermosypnon. 


17.  CONCLUDING  SUMMARY 


As  discussed  in  the  paper  by  Greene,  et  al.  [Ref.  62], 
the  need  for  electric  propulsion  of  minimal  cost,  volume, 
and  weight  is  well  known  and  motivated  by  the  inherent 
advantages  of  this  type  of  propulsion.  A  significant 

constraint  to  be  overcome  in  the  development  and  manufacture 
cf  these  systems  is  the  rate  of  heat  removal  from  these 
devices.  In  1976,  W.  D.  Morris  [Ref.  63]  mentioned  than 
increasing  the  load  per  unit  weight  of  these  devices 
demanded  that  the  inherent  heat  created  by  the  infernal 
losses  be  removed  via  improved  cooling  techniques.  Morris 
also  elaborated  on  the  ’’complex  arrangement"  of  the  ducting 
and  associated  equipment  required  to  accomplish  this 
increased  cooling.  Further,  he  stated  that  "Ther*3  is 
currently  insufficient  technical  information  available  with 
which  to  make  confident  predictions  of  the  effect  of 
Coriolis  and  centripetal  acceleration  on  flow  resistance  and 
heat  transfer  in  the  complex  coolant  circuits...".  Since 
then,  Morris  and  his  co-workers  have  attempted  to  fill  the 
literature  gap  in  this  area,  and  have  uncovered  a  great  many 
facts  regarding  this  flow  structure;  they  have  also  uncov¬ 
ered  a  great  many  questions  for  further  research. 

As  shown  as  early  as  1955  in  an  article  by  Sir  Claude 
Gibb  [Ref.  64],  the  principle  factor  in  the  failure  cf  two 
very  expensive  generating  sets  was  the  inability  of  the 
devices  to  dissipate  their  own  heat.  In  the  late  60's, 
Peach  [Ref.  65]  wrote  of  the  successes  that  the  English 
Electric  Company  had  with  regard  to  liquid  cooling  of  the 
large  generator  rotors  being  used  in  their  500  MW  devices 
and  predicted  that  the  devices  would  reach  1,000  MW  with 
this  technology.  This  was  followed  in  1970  with  a  second 


report  by  Peach  (Ref.  66]  where  he  outlined  the  crrw-h  c  - 
U.S.  generator  sizes  and  predicted  that  the  liquid-coding 
technology  would  be  required  in  the  near  future  *c  support 
the  devices  being  planned  in  this  country  which  exceeded 
1,000  MW.  The  losses  developed  in  the  large-scale,  high- 
power,  supercon duct ive  devices,  currently  being  developed 
for  large-scale  implementation,  require  extensive  use  of 
cooling  at  cyrogenic  temperatures  in  order  to  function  in 
the  super-conductive  state  in  which  they  operate.  A  discus¬ 
sion  of  this  application  is  presented  by  Schwartz  and  Foner 
[Ref.  67],  in  which  a  number  of  the  schemes  discussed  herein 
are  shown. 

The  previous  chapters  have  discussed  the  ability  of 
external,  forced-convective  cooling  to  remove  the 
internally-developed  heat  loads  in  alectrical,  rotating 
machinery.  The  reduction  in  size,  the  increased  efficiency, 
and  the  prolonged  life  cycle  are  all  very  important  reasons 
to  utilize  advanced-cooling  concepts.  The  use  of  turbulent 
flow  is  feasible  at  all  PPM's;  therefore,  both  motor  and 
generator  applications  would  benefit  from  forced-convective 
liquid  cooling  utilizing  turbulent  flow.  Laminar-flow  heat 
transfer  is  feasible  only  at  higher  RPM's,  which  tends  to 
limit  the  application  of  this  method  in  regard  to  motors. 
The  actual  transition  from  laminar  to  turbulent  flow  is  not 
well  known  and  further  investigation  is  required  to  deter¬ 
mine  the  transition  point;  the  transition  may  occur  at  a 
flow  and  RPM  acceptable  for  the  motor  application.  However, 
the  utilization  of  external,  forced-convective  cooling  has 
some  major  difficulties  associated  with  it: 

1)  The  source  of  the  external  cooling  fluid  must  be 
easily  accessible,  dependable  (in  a  casualty-control  sense), 
safe,  volumetrically  small,  light,  and  should  be  reasonably 
inexpensive . 
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2)  The  fluid  itself  should  be  of  high  specific  reef, 
light,  inexpensive,  and  acre  importantly,  ncn-~cxic  and 
non-flammable. 

3)  The  method  of  transporting  the  fluid  from  the  source, 
through  an  external  heat  exchanger,  through  the  device,  and 
back,  should  be  extremely  reliable  and  quite  resistant  to 
external  damage.  T his  is  an  area  of  serious  consideration. 
where  the  rsliablilitv  of  conventional ,  rot  atin.q  sea  Is  is 
questionable  in  the  least. 

The  alternative  to  the  use  of  external,  forced- 
convective  cooling  is  the  utilization  of  a  closed  device: 
(a)  the  closed,  twc-phase  thermosyphon,  (b)  the  closed, 
rota  ting-loop  thermosyphon ,  (c)  the  closed  rotating-loop, 
two-phase  t her mosy phcn,  or  (d)  the  heat  pipe.  These  alter¬ 
natives  remove  the  necessity  of  using  radial  rotating  seals 
in  the  design  of  the  cooling  system. 

The  closed,  two-phase  thermosyphon,  as  shown  in  Figure 
3.6,  must  have  a  sufficient  forced-air  circulation  through 
the  end-bell  regions  to  remove  the  heat.  This  will  require 
optimizing  the'  length  of  the  condenser  sections,  the  geom¬ 
etry  of  the  fins,  the  number  of  fins  and  the  direction  of 
air  flow  within  the  device  to  maximize  the  heat  transfer 
while  minimizing  the  windage  losses.  The  closed  two-phase 
thermosyphon  is  feasible  at  the  higher  HPM's  that  the  gener¬ 
ator  operates  at;  additionally,  this  method  could  also  be 
feasible  in  the  motor  application  if  a  heat  pipe  (i.e.,  a 
thermosyphon  whose  walls  are  lined  with  a  capillary  wicking 
material)  is  used  in  lieu  of  the  two-phase  thermosyphon  for 
the  lower  HPM*s. 

The  closed,  rotating-loop  thermosyphon  requires  an  addi¬ 
tional  heat  •  exchange  through  the  shaft  of  the  device, 
possibly  supplemented  with  forced-air  convection  through  the 
hollow  rotor  region  and  by  the  exposed  radial  sections  of 
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the  loop.  This  shaft  heat  exchanger  may  be  in  the  fcrm  of  a 
rotating  heat-pipe,  as  shown  in  Figures  4.1  and  4.2,  vh-r- 
the  external  exchange  is  to  either  ambient  forced-air,  or  zo 
an  external,  closed  water  spray.  The  shaft  heat  exchanger 
may  also  te  cooled  by  forced  convection  using  liquid 
cooling,  with  axial,  packing-type  seals,  external  to  the 
device  as  shown  in  Figure  4.3.  This  closed-loop  thermcsy- 
phon  could  use  either  single-phase  or  two-phase  operation. 

The  closed  two-phase  th ermosyphon  could  also  be  extended 
along  the  shaft  through  a  bearing  plate  arrangement  and 
cooling  air  sould  te  supplied  external  to  the  casing, 
simplifying  the  internal  configuration,  as  shown  in  Figure 
4.4. 

The  ability  of  these  devices  to  cool  effectively  has 
been  proven  and  their  usage  in  this  application  is  theoreti¬ 
cally  feasible.  T*-e  exact  correlations  for  the  calculation 
of  the  performance  values  for  these  devices,  especially  in  a 
rotating  reference,  will  require  further  research. 

The  analysis  herein  has  shown  the  ability  of  the  coding 
schemes  discussed  to  remove  the  requisite  heat  load  in  order 
to  extend  the  life  cf  the  insulation,  increase  the  effi¬ 
ciency,  decrease  the  weight  and  size  of  motors  and  genera¬ 
tors.  This  study  graphically  presents  the  applicable 
correlations  for  the  determination  of  the  heat  transfer 
within  a  rotating,  electric  device.  It  also  discusses  the 
limitations  of  these  correlations.  It  has  been  shown  that 
closed,  two-phase  systems  of  an  inherently  higher  reli¬ 
ability  are  feasible,  and  deserve  further  evaluation. 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


It 


A.  CONCIOSIONS 

Eased  upon  the  information  reviewed  and  analyzed  in  this 
thesis#  the  following  conclusions  are  made: 

1)  Forced-convective  liquid  cooling  is  a  practical 
means  for  improving  the  efficiency#  reducing  the  size  and 
weight#  and  extending  the  life  of  high-power  motors  and 
generators.  The  main  disadvantage  of  forced-convective 
liquid  ceding  is  the  requirement  of  radial  rotating  seals. 

2)  Closed-loop  liquid  t hermosy phons  are  feasible  at 
high  BPM.  They  offer  improved  reliability  over  forced- 
convective  liquid  ceding  which  utilizes  rotating  seals. 
However,  this  method  requires  the  use  the  shaft  as  a 
secondary  heat  exchanger. 

3)  The  closed,  two-phase  thermosyphon  is  feasible  at 
the  higher  RPMs  of  the  generators.  This  method  requires 
forced-air  convection  through  the  end-bell  regions  or  the 
extension  of  the  devices  along  the  shaft  for  external 
forced-air  convection,  both  of  which  will  add  to  the  windage 
losses. 

4)  Per  low  RPM  applications,  such  as  motors,  the  use  of 
a  heat  pipe  will  be  required  to  overcome  the  inability  of 
the  lew  acceleration  field  to  transport  the  cooling  fluid. 
Both  the  closed,  twe-phase  thermosyphon  and  the  heat-pipe 
require  forced  convection  at  the  condenser  ends  remove 
the  heat  either  internal,  or  external  to  the  device.  This 
may  add  to  the  windace  losses. 
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5)  Correlations  and  theoretical  analyses  are  r.ct  avail¬ 
able  in  the  literature  for: 

a)  the  ccun  ter-curr  snt  entrainment  limit  for  rvn. 
phase  therracsv phcns  under  rotation  and 

b)  the  transition  from  laminar  to  turbulent  flow 
for  rotating  systems. 

E.  BECCMHBNDATIONS 

The  following  recommendations  to  continue  this  work  are 
made : 

1)  Testing  of  the  forced-convective  liquid  cooling 
scheme  should  be  accomplished  using  the  DTNSRDC  test  rig  to 
confirm  the  analysis  and  the  validity  of  the  available 
correlations,  or  to  obtain  proper  correlations  for  the 
specific  geometry  and  utilization  being  projected. 

2)  An  experimental  program  should  be  devised  and 

conducted  to  determine  the  entrainment  limitation  in  a  hori¬ 
zontal,  heated,  channel  rotated  about  a  parallel  axis  with 
counter-current  flow  cf  liquid  and  vapor. 

3)  An  experimental  program  should  be  developed  to 

determine  the  optimum  condenser  length  and  proper  fin 
arrangment  for  the  use  of  both  heat  pipe  and  closed  two- 
phase  thermosy phcn  systems  for  rotor  cooling.  Consideraf ion 
cf  the  structural  integrity  must  be  included  for  the 
condenser  end-lengths. 

4)  An  experimental  program  should  be  developed  for 

correlating  the  laminar- to-turbulent  transition  during 
forced-convective  liquid  cooling  in  a  heated,  horizontal, 
rotating  pipe. 

5)  An  experimental  program  should  be  developed  to 

modify  the  inside  geometry  of  the  closed-loop  thermosyphon 
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(and  axial,  closed,  two-phase  thermcsyphor.)  to  optimiz-  th- 
heat-transf er  and  the  friction  losses  wi-hin  the  device. 

6)  A  detailed  analysis  should  be  made  of  the  shaft 
cooling  potential  for  the  rotating- loo p  thermosyphon  with 
regard  to  the  heat-transfer  enhancement  cf  the  shaft  to  the 
fluid  and  the  mechanical  design  of  the  shaft. 

7)  An  experimental  program  should  be  developed  to 
analyze  the  closed,  two- phase  thermosyphon  with  radial 
condenser  sections  and  the  closed- loop  two-phase  thermosy¬ 
phon.  These  two  devices  may  experience  start-up  balance 
problems. 

8)  Various  configurations  should  be  analyzed  for  their 
potential  for  heat  transfer  to  include  a  combination  of  the 
systems  described  herein,  as  depicted  in  Figures  4.1  through 

4.4. 

9)  The  concept  cf  modifying  the  geometry  of  the  rotor 
windings  within  these  devices  in  order  to  minimize  heat 
conduction  path  lengths  to  the  casing,  shaft,  etc.,  and 
thereby  the  internal  temperature  gradient,  should  also  be 
investigated. 


APPENDIX  A 

GOVERNING  EQUATIONS  FOR  LAMIHAR  CONVECTION  IN  UNIFORMLY 
SEATED  BORIZONTAL  PIPES  AT  LOB  RAYLEIGH  NUMBERS 

The  following  development,  by  Morton  [Ref.  11]*  i^ 
presented  herein  for  completeness.  The  steady  laminar 
moticn  cf  a  fluid  in  a  horizontal  circular  pips  of  radius  a, 
the  walls  of  which  are  heated  uniformly  so  that  a  constant 
temperature  gradient  t  is  maintained  in  the  direction  of 
the  axis.  The  flow  will  be  referred  to  in  cylindrical  coor¬ 
dinates  (r,<f>  ,  z)  with  4>  measured  from  the  upward  vertical 
and  the  z-axis  along  the  axis  of  the  pipe;  the  velocity 
components  are  denoted  by  u,v,  and  w.  The  effects  of  dissi¬ 
pation  and  cf  the  pressure  term  in  the  energy  equation  will 
be  neglected  and  that  variations  in  the  density  due  to  the 
temperature  differences  are  so  small  that  they  only  effect 
the  buoyancy  term.  Thermal  conductivity  and  kinematic 
viscosity  are  assumed  to  be  constant  (which  will  introduce 
quantitative  errors  into  the  solution,  but  should  not  change 
its  general  character) . 

The  equation  of  continuity  is  then: 

3ru/3r  ♦  3v/3<j)  ♦  3rw/3z  =0,  (A - 1 ) 

and  the  energy  equation  is; 

u  (3T/3r)  +  v/r(3T/3<|>)  ♦  w  (3  T/3  z)  =  k  V^T,  (A. 2) 

where  V*=3*/3r«  +1/r(3/3r)  ♦  1/r*  (3  */3<j>  *)  +3V3z2  and  T  is 

the  local  fluid  temperature.  The  momentum  equations  can  be 
written  in  the  form  cf  equations  (A. 3)  through  (A. 5). 


98 


u  (  dU/  a!)  +  v/r  (  3u/  3  i)  +  w  (  3u/  3z) -v  */z  =  -1/p(5?/5£) 


v(  V2u-u/r2 -2/r  2  <  3v/ 3i)  -  a  (T  -T)cos;>  ,  (A. 3) 

w 

u  (3v/3r)  +  v/r  (  3v/ 3$)  +w (  3v/  3z)  =uv/r  =  -  1/pr  (3p/3>J>) 

7  (  72v  +  2/r2  (3U/3  J-v/r2)  +  g(Tw-T)sir£  ,  (A. 4) 

u  (  3w/3r)  v/r  (  3w/3$) +w  (  3w/ 3z)  = 

-1/p  (3p/3z) +v  v2w*  (A. 5) 

In  equations  (A.  3)  and  (A .  4)  ,  the  buoyancy  force  has  been 
calculated  relative  to  the  fluid  at  the  same  level  adjacent 
to  the  pipe  wall  and  the  remaining  distribution  of  force  has 
been  absorbed  into  the  pressure,  p. 

For  steady  convection  sufficiently  far  from  the  pipe 
opening  to  avoid  inlet-length  effects,  the  temperature 
throughou*  the  pipe  increases  uniformly  with  the  distance 
along  the  pipe  axis.  Hence,  the  distribution  of  buoyancy 
force  in  sections  of  the  pipe  is  independent  of  z;  as  the 
secondary  flow  is  caused  by  the  buoyancy  forces,  the  flow 
field  must  also  be  independent  of  z.  It  follows  that  there 
should  be  a  similarity  solution  with  u,v,w,  and  TW~T  as 
functions  of  r  and  0  only,  and  with  p=Yz+P(r,  4>  )  (where  P 
contains  the  terms  absorbed  into  p)  . 

The  continuity  equation  reduces  to  the  form: 

3(ru)  /3r  +  3v/3<j>  =  0; 

hence,  a  dimensionless  Stokes  stream  function  <p  can  be 
introduced  in  such  a  way  that: 

ru/v  =3iJj  / 3 4>,  v/v  =  -3ip/3r. 
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The  main  aquations  (A.  1)  through  (A.  5)  may  be  racaroi  no 

non -dimensional  form  by  f ha  transformations  r=aR,  y=(v/a.13, 

and  T  -T=taPr  ,  where  Pr  is  the  Prandtl  numbar  (v/cc)  .  If 
w 

the  pressure  is  eliminated  between  equations  (A.  3)  an 3 
(A.  h )  ,  the  resulting  momentum  equations  are: 

7*  t  ♦  V?  (3t/asO/3y)-ar;-/?oO/?R)  }  VZ'P  * 

Ba{(39/3H)  sin®  ♦  1/R  (30/3 i) cos  j>},  (A. 6) 

VfW+  1/Bf(3T|;/3R)  (3/3(J>)- 

(dip/d<p)  (3/3  R)  }  W  +4Re  =0,  (A. 7) 

and  where, 

Vf  =  3  2/3R 2  ♦  1/R  (3/3  R)  +  1/R2(  3  2/3cj>  2)  . 

The  Rayleigh  number.  Pa,  is  given  by  gBia*/(av  )  and  the 
Reynolds  number  is  the  normal  Reynolds  number.  Re,  for 
Poiseuille  pipe  flow  based  on  the  pipe  diameter  and  the  mean 
velocity  across  a  pipe  section  is  (a/v)  (  ya2/4y  )  . 

The  solution  of  the  energy  equation  depends  on  the  temp¬ 
erature  cn  the  boundary  condition  at  the  wall  of  the  pipe. 
Pipes  will  usually  have  reasonably-thick  walls  of  material 
with  thermal  conductivity  much  higher  than  that  of  the 
fluid,  so  there  will  be  little  variation  in  wall  tempera¬ 
ture  around  the  pipe  circumference.  This  will  be  specially 
so  for  slow  rates  of  heating  when  asymmetries  of  the  flow 
will  have  small  amplitude.  Hence,  write  T  =T0+  z,  where  TC 
is  the  wall  temperature  in  the  section  containing  the 
origin.  Although  there  is  uniform  heat  transfer  per  unit 
length  of  pipe,  the  local  heat  transfer  will  be  slightly 
greater  near  the  bottom  of  the  pipe  than  near  the  top.  with 
this  assumption,  eguation  (A. 2)  reduces  to  the  non- 
diaensional  form  of  eguation  (A. 8). 
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72  a  ♦  (a/E)  {  (  9^/3r)  (3  9/ 3(D)  -(Z\p/d<p)  (36/3r)  }  +w=0 

The  boundary  conditions  are: 

u,v,v,  and  3  are  zero  on  r=a; 
u,v,wr  and  0  are  finite  on  r=0. 

(A. 9) 

Although  solution  of  equations  (A. 6)  through  (A.  3),  satis¬ 
fying  conditions  of  (A. 9  is  considerably  difficult,  succes¬ 
sive  approximate  3  to  the  solution  can  be  obtained  by 
expanding  'l>,  W,  nd  9  as  power  series  in  the  Rayleigh 
number,  provided  that  this  is  numerically  small.  Supposing 
that : 

ib=Ra^.+Sa2^  +..., 

1  2 

w=w  0+  BaW  !+Ra  2w  2+. . . , 

0= 6 0 ♦ Ea 9 j +Ra 2 9 2 + . . ,  . 

(A. 10) 

The  leading  term  ’P0  cf  ip  must  vanish  because  there  is  no 

circulation  when  A  is  zero,  but  60  is  not  zero  since  the 

difference  in  temperature  between  a  fluid  element  and  the 

neighboring  wall  is  proportional  to  x9  .  When  relations 

(A.  10)  are  substituted  ir.  equations  (A.  6),  (A. 7)  ,  and  (A.  8) , 

three  sets  of  eauations  for  the  funtions  ip  .  w  ,  and  9. 

i  i  1 

are  obtained  by  equating  coefficients  of  powers  of  Ra. 

From  equation  (A. 7)  ,  the  basic  equation  is: 

V2K0  ♦  URe  =  0,  (A.  11) 

which  has  a  solution: 

W  Q=  Re ( 1-r2)  (A.  12) 

which  satisfies  the  conditions  W  0=0  at  5  =  1,  w0  finite  at 
R=0.  This  is  ordinary  Poiseulle  flow  under  a  pressure 
gradient  -4  p^Re/a3  in  an  unheated  pipe. 


The  first  of  the  equations  derived  from  equation  (A. 8) 
for  the  temperature  distribution  is: 

29  +w  =o,  (A-  13) 

V  x  v  0  ' 

subject  to  9  O=0  at  5  =  1,  9,  finite  at  R  =  0,  and  the  solution 

is  axially  symmetrical  as  no  account  has  been  taker,  of 
gravity  at  this  stage  of  the  approximation.  This  is  satis¬ 
fied  by: 

0O  *0/16)  (3-R*)  f  (A-  14) 

which  is  the  customary  solution  for  forced  convection. 

The  first-order  approximation  for  \p  satisfies  the  equa¬ 
tion: 

7*-^  =  (390/33)  sin  ■$  (A.  15) 

obtained  as  the  coefficient  of  Ba  from  equation  A.  6,  and  the 
boundary  conditions  3’4>i/34>=0  on  R  =  1 ,  and  R~l  (3lN  /«*<!>)  , 

3(j,i/3  r  remain  finite  at  R=0.  If  'Pl  is  assumed  to  have  the 
form  i|;  j  (B)  sin<+>  ,  the  dependence  of  equation  A. 15  on  <j>  is 
eliminated,  and  ij>  (H)  is  easily  found: 

^*-(1/4608)  ReR(1-E2)  2(10-R  «)  sin  <P  .  (A.  16) 

Similarly,  the  rfrst-order  approximation  to  (^satisfies  the 
equation : 

72(1  =  (1/P)  (3j»  j/94>)  (9W0/3R)  ,  (A. 17) 

i  i 

and  the  boundary  conditions,  W : -0  at  R=1  and  finite  at 
R=0.  The  dependence  of  equation  A.  17  on  <i>  is  eliminated  by 
faking  M^HjOUcoso  ,  whence 


W  .  =  -  (1/1  3  432  0)  Re2R  (1-R2)  (49-5  1R2«.  iqr  +  -*6)  co s4>. 


(A.  1Q) 


The  f ir st- or der  approx i car ion  to  9  satisfies  th-  aqua¬ 


tion: 


V20  ,  =  (o/R)  (3<P,/3<fr)  (30  /3  R)  -H  , 

j  1  1  1 


(A.  19) 


and  the  boundary  conditions,  ?  t  =9  at  r-1  and  :  finite 

B=0.  Hence, 

0  !=-  (1/22118400)  Re  2g(l -32)  {(38  1  + 1325a )  -(354 +1675a)  R2  + 


+  (146+9250  )  R*-(29  +  200c)  R«+(1-2a)  K8}  cos©  . 


(A.  20) 


This  completes  the  approximate  solution  to  the  first 
order  in  Ra  for  convection  in  a  horizontal  pipe,  which  is 
heated  so  that  there  is  a  constant  temperature  gradient 
along  the  walls  and  uniform  temperature  around  the  girth. 
It  may  he  noted  frcm  equations  (A.  10),  (A. 12),  (A. 14), 

(A.  16),  (A.  18),  and  (A.  20)  that  the  full  convection  solution 
depends  essentially  cn  the  product  RaRe.  This  is  clear  if 
it  is  recalled  that  Ra  is  proportional  to  the  increase  in 
wall  temperature  along  the  pipe  and  Re  is  proportional  to  a 
characteristic  velocity  of  the  flow  along  the  pipe.  An 
increase  in  Ra  means  that  the  fluid  will  be  carried  through 
a  larger  temperature  differences,  and,  hence,  there  will  be 
increased  buoyancy  forces;  but  the  same  effect  can  be 
produced  by  an  increase  in  Re,  and  hence  in  the  velocity  of 
the  main  flow. 

The  two  most  interesting  consequences  of  such  a  solution 
are  the  modification  of  the  heat  transfer  due  to  secondary 
flow,  and  the  effect  of  heating  on  the  tangential  stress  at 
the  wall  (or  on  the  rate  of  flow  for  a  given  pressure 
gradient).  However,  as  both  w  ,  and  9,  vary  with  cost))  , 


neither  the-  flux  ncr  the  heat  transfer  across  the  v.cl  •. 
section  c f  the  pips  will  fce  chnncad  by  the  firs^-tr ’-r 
aporeximation  a  1  •‘■houah  there  will  be  local  variations. 

Tc  find  *he  changes  ir.  total  flux  and  heat  transfer,  it 
is  necessary  to  proceed  to  the  second-order  approximations, 
and  the  revelant  equations  for  these  are: 

r*.j,  =  (1/R)  7?*.  ,<!>■) /3(R,d>  )} 

i  -  *  i  i 

♦(36/3  3)  sin**  (1/R)  OQ/34>)coS(j,  , 

=  (1/a)  0(-W  ,'l>  )/3  (a,  o)}+  (1/H)  [3(»  ;p  J/3(R,  4)}  , 

12  11  0  z 

72 a  =(a/R){3(0  ,ip  ) /3  (B,$)  }  + (a/R)  (3  (6  , ip  /3(R,$)}-W  . 

ju2  02  2 

(A. 21) 

The  only  difficulty  in  solving  these  equations  is  numerical 
tediom  and  the  remaining  second-order  solutions 
presented  in  [Ref.  11]. 
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&E21SDIX  3 

SINGLE-EHAS2  ANALYSIS  (COMPUTER  PROGRAM  AND  RESULTS) 

The  following  program  was  written  to  give  the  results  o 
the  various  correlations  for  the  heat  transfer  for  th 
force!,  single -phase  convection  in  the  rotating  reference  a 
disscussed  in  Chapter  II.  The  program  is  so  written  tha 
modification  to  include  other  correlations  is  quite  easy 
It  was  written  in  HP-Basic  for  the  HP-9826  computer 
Sub-programs  for  the  graphics  are  not  shown.  Sample  run 
are  included  immediately  following  the  listing. 

-Pages  105-109  Program  Listing 

-Pages  110-111  Thermophysical  Properties  (Common  to 

Appendices  3  through  D) . 

Sample  Calculation 


-Page  112 


1  Z.J  • 

vd.riu.ai  y  _  .  E1 0  *■’ 

_ru*  •  -\r  '* .  t) : 

■ '  4  0 

COM  'Cm  in/  Xmin . Tram . Sf x .Sf y 

i  050 

optMTr^  * 

060 

BEEP 

•070 

pppr  USING  " 

2X  ,  De*  aij.lt  values:  . 

1080 

D*4 . 763E-3 

!  Tune  aiamerer  (n> 

090 

j- . 8230 

1  ' <jDe  length  (m> 

HOC 

R- .4021 

1  R aa  i  as  o  f  rotor  ( n ) 

1110 

PRIM'  USING  1 

4X,  :  one  diameter  -  "".Z.4DE,""  <m> . :  D 

1120 

PRINT  USING  1 

4X  .  ""Tube  length  =  ""  .  DD.D  <m> . 

; '  3C 

PRINT  USING 

4X . ""Radius  of  rotor  -  "".Z.  3D,""  in) . :.R 

1140 

BEEP 

■50 

INPUT  "OK  TC 

ACCEPT  DEFAULT  VALUES  < 1  - Y . 0-N > 7" . Id 

1160 

BEED 

i ;  70 

INPUT  "LIKE  A 

h«RD  COPT  ('»Y,C*N)7,,.Ihc 

1180 

IF  Ihc*1  THEN 

1  ’90 

PRINTER  IS  70 

* 

1200 

E:_$E 

i  2’  0 

°R INTER  IS  1 

1220 

END  IP 

1230 

IF  Id-1  THEN 

1300 

1240 

BEEP 

1250 

INPUT  "ENTER 

TUBE  DIAMETER  (m>".D 

1260 

BEEP 

'■  270 

INPUT  "ENTER 

tube  length  <m)".L 

1280 

BEEC 

•  ‘■40 

INPUT  “ENTER 

RADIUS  OF  ROTOR  <m>".R 

1  300 

DRINT  USING  ‘ 

1 0X. Geometric  Variables  —*  """ 

•3:0 

PPJNT  USING  1 

’4X.""Tu&e  diameter  <D>  *  "" . Z . 3DE <n) . :D 

•320 

PRINT  USING  1 

14X .""Tube  length  <l>  -  "" . DD  .  DD  .  ""  <  n > . :  l 

1  130 

PRINT  USING  " 

14X.""L/D  -  ""  . 5D  . D"  : •_ / D 

1340 

PRINT  USING  ’ 

14X. ""Rotor  radius  <R)  -  . Z  .  30  .  (m) . :R 

•  350 

PRINT 

•360 

PRINT  USING  • 

10X.  .omputed  values  are: 

’370 

printer  IS  1 

’380 

BEEP 

1390 

PRINT  USING  1 

2X. ""SELECT  OPTION: . 

1400 

PRINT  USING  ‘ 

4X .  1  Single  point"" 

1410 

PRINT  USING  1 

4X.""2  Nu  versus  RPM . 

M20 

PRINT  USING  ’ 

4X .  3  Nu  versus  Re""' 

1430 

PRINT  USING  1 

4X .  4  Nu  versus  Prod'”1 

1440 

INPUT  Qp 

1450 

IF  Op >2  THEN 

1460 

PRINT  USING  1 

4X  , ""SELECT  OPTION: . 

1470 

PRINT  USING  1 

4X.""2  VARY  RPM . 

1480 

PRINT  USING  1 

4X,"" 3  VARY  Re . 

1490 

INPUT  Qpx 

'500  IF  Gox-2  tHEN  Op-2 

1510  IF  Opx-3  THEN  Gp-3 

1520  END  IF 

1530  IF  Op>1  THEN 

1540  BEEP 

1550  J’*0 

•560  PRINTER  IS  ' 

’570  °RINT  USING  "4X . ""SELEC'  OPTION: . 

1580  PRINT  USING  "4X.’”‘  1 -Di t tus-&oe : ter .  2-Nakayama,  3-Nakavama 


vSiNO  "4X, 

r  *  H  (if-  <d  v  .iPt  oL  “  - 

INPUT  Ian 


J-iooas  a  “orr:-:  .  jt?cnenson , 


iiiruT  '_In£  7‘J  r.ui  i  1*1  .u*N)  '  '.Jkolot 
IF  Okpiot-1  THEN 
CALL  P lot ( J j . Jc . jd . X . Y . Type .Cx > 

3  EEP 

INPUT  "SOLID*? .DASH-DASH-2. DOT-DOT. 3". Type 
END  IF 
END  IF 
Nstep* 1 00 
IF  0o-2  THEN 
BEEP 

IMP'JT  "ENTER  RPM  RANGE  U1IN.MAX)’’  .Rpml  .RomH 

Rpm-Rpm 1 

ELSE 

BEEP 

INPUT  "ENTER  RPM  OF  ROTOR". Rpm 

END  IF 

BEEP 

INPUT  "ENTER  OPTION  < ? -0 ,2-Mf . 3-Dt )".Io 
IF  Ihc-1  THEN  PRINTER  IS  701 
IF  Ihc-0  THEN  PRINTER  IS  1 
PRINT 

PRINT  USING  "1  OX, '”’***  Operating  Variables  *** 
PRINT 

PRINT  USING  " I  OX . ""Input  variables  are:  . 

IF  Io-1  THEN 
BEEP 

INPUT  "ENTER  COOLANT  NASS  FLOU  RATE  <lcg/s>"  .Mf 


DRINT  USING  "  1  AX  . '"'Coolant  mass  rate  *  ’"' .Z.3DE."" 

INPUT  "ENTER  INLET  AND  OUTLET  TEMPS  <C>’’.Ti.Te 
PRINT  USING  "? AX . ""Coolant  inlet  temp  «  "“.DD.DD."" 

PRINT  USING  "1 4X . ""Coo lant  outlet  temp  =  ’"’.DD.DD,”" 

END  IF 

Ic  Io-2  THEN 
BEEP 

INPUT  "ENTER  HEAT  LOAD  <U)".Q 

PRINT  USING  ”1  AX  .’"'Heat  load  *  .  4D.DD.'" 

BEEP 

INPUT  "ENTER  INLET  AND  OUTLET  TEMPS  <0".Ti.To 
PRINT  USING  ” 1 4X , ””Coo i an t  iniet  temp  *  '"'.DD.DD.’” 

PRINT  USING  "? 4X , ""Coolant  outlet  temp  *  '"’.DD.DD."" 

END  IF 

Ir  Io*3  THEN 
IF  Jj>0  THEN  2200 
BEEP 

INPUT  "ENTER  HEAT  LOAD  <W)’’.0 

PRINT  USING  ”14X  ,  ’"’Heat  load  *  '"’ .  4D  .  DD  ,  ’" 

BEEP 

IF  0p<3  THEN 

INPUT  "ENTER  COOLANT  MASS  FLOW  RATE  <kg/s>".Mf 

PRINT  USING  "1 4X . '"'Coolant  mass  f lou  rate  -  "".Z.4DE,"' 

ELSE 

INPUT  "ENTER  MASS  rL0U  RATES  < MIN  .MAX)" .Mf 1 . Mf h 
Mt-Mf  1 
END  IF 
BEEP 

INPUT  "ENTER  COOLANT  INLET  TEMP  (C>".Ti 


( if  g/s )  ’ 


<c> . 

<c> . 


. 4D.DD."" 

"".DD.DD.*"’ 

•"’.DD.DD."" 


’’’ ,  4D .  DD  , 


(k  g/s  > ' 


22' 0  END  IF 

Z<-<. .1  V>  iu*.  u'.j 

2230  Cp”I-  NLpu (  Ta> 

2240  Mu*FNt1uw(Ta> 

2250  K«F NKw(Ta) 

2250  sr*cNPru<Ta) 

2270  Rho“FNRrtou(  la.) 

2230  Nuv-Mu/Rho 
■'"'•41  rr  t,  =  .  TH^N 

2300  Q-Hf«CD*(To-Ti> 

2310  PRINT  USING  "  1 4X  , ""Heat  load  «Q>  *  .4D . D  . <U> . :0 

2320  END  IF 

2230  IF  Io-2  THEN 

2340  Mt-0/(Co«(To-Ti>> 

2250  PRINT  USING  "! 4X . ""Coolant  mass  flow  rate  -  "".Z.4DE.""  <kg/s> . :Hf 

2360  END  IF 

2370  IF  Io-3  THEN 

2380  Toc”Ti+Q/<Mf«Cp> 

2390  IF  PBS<  To-Toc)>  .01  THEN 
2400  To* < To+Toc )» . 5 
2*10  GOTO  2220 
2420  END  IF 

2430  IF  Jj*0  THEN  PRINT  USING  "MX. ""Coolant  outlet  temp  -  . 3D  . DD  . ""  <C> 

“ :  To 


2440  END  IF 
2 450  IF  J?-0  THEN 
2460  PRINT 


2470 

PRINT 

USING 

•M0X. 

""Fluid  properties  eva 

1 uated 

at  "".DD. 

DD.  ""  <C> 

are  :  M* 

'  d 

2  4  9  0 

PRINT 

USING 

”1  4X . 

""Soecific  heat  <Cp> 

*  ' 

. 4D . D . ""  ' 

< J/kq.K) . 

:Cp 

2490 

PRINT 

USING 

"MX. 

""Viscosity  (Nu> 

K 

. Z . 4DE  . "" 

<  N  .  s  /  m  ’  2  > " 

;Mu 

2500 

PRINT 

USING 

"14X. 

""Tnermal  cond  <k  > 

=  ' 

. Z . 4D . 

<W/n.K> . : 

K 

2510 

PRINT 

USING 

"14X. 

""Prandti  number  (Pr> 

* 

. Z . 3D" : Pr 

2520 

PRINT 

USING 

"14X. 

""Density  <Rho> 

X 

.40.0."" 

(kg/m  3> . 

:  Rho 

2530 

PRINT 

USING 

"14X. 

""Kinematic  vis  <Nuv> 

=  •••* 

. Z.4DE."" 

( m  2/s  1 . 

:  N  uv 

2540 

9eta*FNBeta<  Ta) 

2550 

PRINT 

USING 

"  1  4X  . 

""Coef  ther  exp  <Beta> 

_  i»»i 

. Z . 4DE . "" 

(  1 /K )  : Beta 

2560 

PRINT 

2570 

PRINT 

USING 

"  1  OX . 

'"'»«*  Calculations  *••** 

2580 

PRINT 

2590 

PRINT 

USING 

"  1 0  X . 

'"'Prel  immary  calculations:"" 

2600  END  IF 

2510  Re-4»Mf/(PI*D*Mu> 

2620 ! 

2530!  Friction  Factors  for  stationary  reference 
2640  IF  Re<2*1 .E+4  THEN  2670 
2650  F- . !84/Re  '  .2  !!4DEqn8.21 

2660  GOTO  2630 

2670  F* .316/Re' .25  ! I4D  Eqn  8.20 

2580! 

2530  0pp-Q/(PI*D«L) 

2700  0p*0pp*PI*D 

27  T  0  Vm-4*Mf/(Rho*PI*D'2) 

2720  Dps*F*Rho*Vm"2*L/<2*D) 

2730  Gmega*Rom*2»PI /60 
2740  Ro*Vm/<Omega*D> 

2750  Jay*0mega*D '2/Nuv  !From  STEPH-Rotat lonal  Re 
2760  Kw*FNKw( Ta> 

2770  Ot t*  <  To-T  i  > /L 


2730 

2790 

2300 

CO  •  vj 

232U! 

2830! 

2840 

2850! 

2S60 ! 

2870 

2380! 

2890! 

2900 

er 

2910! 
2920! 
2930 
blent) 
2940 
2950! 
2980! 
2970 
2980 
2990 
3000 
3010! 
2020 ! 
3030 
3040 
3050! 
2060! 
3070 
3080 
3090! 
3100! 
3110 
3 ’20 
3130! 
i  140! 
3150 
3160 
3 '70 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
2250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 


Gr-R*<0mega  2>*3eta*<  <D/2>  4)«Dtt/Nuv  2 
Ra-Gr*pr 

6amroa-<Re  < 22/1 3> ) •< ( Gr- (Pr  .b>>  <-12/!3>> 

rv  t  *  n  C  \ 

1  y “  1  /  v  _  .  w  ' 


Laminar  Value 
Nulam*48/1 1 

Turbulent  Value  { rom  Di t tus-3oe 1  ter 
Nudb* . G23»Re~ . 8*Pr '  .4 

Stepnenson  Correlation 

Nus- < P r 4/ . 367 ) » . 007 1  -Re ' . 88* Jay  . 023  ! Ai r . Turb_E3N2 .23-correc ted  tor  uat 


Nakayama  Correlation  for  '.urDulent  condition; 

Nunak-<Re' .8>*<Pr‘.4)*<Tg <1 /30>  >-H  - .  0"  4/<  Tg>  (1/6>>-.033  !  Eqn  2.9  <Tur 

Ty-Nunak / (Re ' . 8*Pr ‘ .4) 

Nakayama  Correlation  for  iaminar  conditions 
J2- J-2 

Sq-2/1 1*< 1 +SQR< 1+77/4«< !/Pr'2>> ) 

Cf«1-.486«<  <3«Sq-1 >‘ .4/<Sq-<Sq-Pr-SQR(5>+2>>  >-J2/(Ra-Re>'  .6 
Nunak 1-48/ 1 1  * . 1 91 /Sa*<3*Sq- 1 > ' .2»<Ra*Re>  '  .2*’ /< ’ +( ! /10-Sq-Pr>  > 

Nakayama/Fu: toka  Correlation  for  Radial  Pioes 
Nunf-<  .014/.023)*<Re/Ro'2.5)',.124*Nudb  !  TURB-EQN  2 .21 

Mus-Mu  !Temp  Value  Mus'jMu 

Sieder-Tate  Correlation  for  Turbulent 
Nus  t- . 027*Re  8»Pr ' . 333 3»( Mu/Hus >.14 
Prod-Ra*Re-Pr 

Woods-Horr is  Correlation  for  Laminar 
Nuum* . 262*Prod' . 1 73-48/ 1 1  1  Eon  2.22~Wtt2 
Jl -Jay/8 

Woods-Morr is  Correlation  for  Radial  pipes 
Nuwm2- . 01 5-Re~ . 78-J1 '  .25  ! f rom  um2 
IF  Jj-0  THEN 

pRINT  USING  " 1 4X , '"’Reyno Ids  nurnDer  <Re> 

PRINT  USING  "14X  .  ”**Fr  ict  ion  factor  <stat> 

PRINT  USING  ” 1 4X , ""Heat  flux  <Qpp> 

PRINT  USING  "1 4X . ""Mean  fluid  *el  <Vm> 

PRINT  USING  ”14X.”"RPH 
END  IF 

IF  Ihc-0  THEN  PRINTER  IS  1 
PRINT 

IF  Jj-0  THEN  PRINT  USING  ’’  1 0X  . ’’"Resul  ts  : 

IF  Okplot-O  OR  <0kplot-1  AND  Nstep<11>  THEN 

PRINT  USING  “10X .  Nudb  Nunak  Nunf  Nuwm  Nus  Nuuim2 

PRINT  USING  “1 3X , 6 (3D . DD . 2X )” :Nudb . Nunak . Nunf . Nuwm . Nus . Nuum2 
PRINT  USING  ”  1 4X , Ra  Ro  Ra-Re-Pr  (if  Nak  1  ” 

PRINT  USING  ”14X.4(Z.2DE.2X) .3D . DD ,DD . DDD . Z . DD” :Ra. Ro .Prod . Mf . Nunak 1 

PRINT  USING  " 1 4X  , RaRe  Tg  Ty . 

PRINT  USING  ”13X , D . 2DE .2X . Z . 3DE .2X . Z . DDD" :Ra»Re . Tg . Ty 
END  IF 

IF  Okelot-1  THEN 
IF  Ian-1  THEN  Y-Nudb 
IF  Ian-2  THEN  Y-Nunak 


-  .  Z .  3DE"  :Re 

-  "".Z.3DE”;F 

-  , Z .  3DE (W/m  '2) . :Qpp 

-  ""  ,  Z . 2DE . (m/s) . :  Vm 

-  "" . 40 .DD : Rpm 
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'  C*  1  0  DIH  k ( 8  > 

•  1'.'.  "'£-£)  -  T  ••-■i'1  ^co:  _  •  ;  t  • -''.cc  JiC  c  >  . '%32c  c'"/l  .1 

•  C  ju  Dn  m  m  .  ■  o71  i  73i  .  2C  .  573vj(  .  :■ 

'040  READ  K  <  *  > 

1  050  T- < Tsteam+273 . 1 5> /647 . 3 

TOGO  Suin*0 

1070  FOR  N-0  TO  4 

1  080  Sum”Sum+K < N > * < 1  - T  )  (N+1) 

1090  NEXT  N 

’  ’  00  Br'Sum/f  T*<  1+K<5>*(!-T>+<(G>*(,-t)'2))-<’-T>/(K'^7)*(1-T)2+i<<8)> 
'I'O  Pr*EXP(3r ) 

1 1 20  P*22 1 2000G*Pr 

1130  RETURN  P 
H40  FNEND 
1150  DEF  FNHf  g<  T  > 

1 'GO  Hfq«2477200-2450*(T-i0> 

1170  RETURN  Ufa 
1180  FNEND 
1190  DEF  FNMuu<  T  > 

1 200  A-247 .8/<T+133. 15) 

1210  Mu-2.4E-5*10 "A 

1220  RETURN  Mu 
1230  FNEND 
1240  DEF  PNVvst<Tt> 

1250  P-FNPvst ( T  t ) 

1260  T-T  t+273 . 1 5 

1270  X-1500/T 

1280  F1*1/< 1+T*1 .£-4) 

1290  F2-  <  1 -EXP(-X > ) ~2 . 5»EXP< X ) /X'  .5 

1300  B- . 00 1 5*F 1  - . 000942»F2- . 0004882*X 
’310  K«2*P/< 461 . 52»T ) 

'  320  .5>/K 

1330  RETURN  V 

•  340  FNEND 

1350  DEF  FNCoulT) 

1 360  Cdu-4.21 120858-T«(2.26326E-3-T*<4.42361E-5*2.7l428E-7*T>) 

1370  RETURN  Cowl  000 
’ 380  FNEND 
1 390  DEF  PNRhoui(T) 

140C  Ro * 999 .52946 +  T*<  .  0 1 269-T-*  <  5  . 4825 I  3E-3-T-*  1  ,234147E-5>  > 

1410  RETURN  Ro 

1420  FNEND 

1430  DEI7  FNPrui<T> 

i 440  Prw»FNCpw<  T>*FNMuw< T ) /FNKu< T ) 

1450  RETURN  Pru 
1460  FNEND 
1470  DEF  FNKu(T) 

1480  X-1T+273. 15>/273. 15 

1490  Kw-  . 92247+X»<2 . 8395 -X*<  1  ,8007-X*<  . 52577- .  07344-X >> > 

1500  RETURN  Ku 

1510  FNEND 

1520  DEF  FNTanh(X) 

1530  P-EXP(X) 

1540  O-EXP(-X) 

1550  Tanh-<P-Q)/<P+Q) 

1560  RETURN  Tanh 
1570  PNEND 
1580  DEF  FNHf(T) 
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'590  Mf-T.<a.2038$9-T«<5.*S!32E-<»-~-a.55150 

1600  RETURN  Hf-1000 

’610  -WEND 

1K20  DE"  FNT./«r(Pi 

'  Oju  ■  U.-  :  1  0 

1690  T1-1C 

1650  Ta*<  Ta+Tl  >■*  .5 

1660  PcmFNPvs t < Ta> 

1670  F  «BS<<P-Pc)/P)>.001  THEN 

1680  IF  Pc<P  THEN  Tl-Ta 

1690  IF  Pc>P  THEN  Tu-Ta 

1700  GOTO  1650 

1710  END  IF 

1720  RETURN  Ta 

1 730  FNEND 

1790  DEF  FNBeta(T) 

1750  Rop*FNRhou< T+  .  1  ) 

1760  Rom*FNRhou< T- .  i ) 

1770  Beta"-2/(Rop+Rom)»<Rop-Rom  > /  .2 
1780  RETURN  Beta 
: 790  FNEND 
1800  DEF  FNAlpha(T) 

1810  Alpha“rNKw<  T )  /<  FNRho'j<  T  >»FNCow(  T  > ) 

1820  RETURN  Alpha 

1330  FNEND 

1890  DEF  FNSigma(T) 

1350  T  t-697 . 3-T 

I860  A* . 00 1 *Tt  *2*< . 1 160936807 /< 1  +  .S3«Tt>) 

1870  B-. 001 121 909688-5. 7528051 8E-6*Tt 

'380  C«1 .28627965£-8-Tt  2-1 . 1 9971929E- 1 1 *Tt 

1890  Siama-A+B+C 

1300  RETURN  Sigma  - 

1910  FNEND 
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•  leone* r  it  harianies  *  *  * 

lore  cianeter  (C1)  *  4 .  7S3E-03  (n) 

**  ,  l_  -  ’ _ !  '  .  ,  _  "'''/-X 

1 -/LI  -  \  ( C  .O 

Rotor  racius  <R>  *  0.402  <m> 

Conoiitec  values  a. re: 

*'*  .berating  Varianies 

Incut  variables  are: 

-feat  1  oac  *  50.00  <w> 

Coolant  inlet  temp  -  05.00  (C> 

Coolant  outlet  temc  *  05.72  <C' 

Coolant  mass  flow  rate  »  1 -S329E-02  (kg/s> 

Fluid  properties  evaluated  at  45,36  <C>  are: 

Specific  neat  (Cp>  *  4224.7  (J/kg.'O 

Viscosity  <Mu)  *  5.8661E-04  <N.s/m‘2> 

Thermal  cond  <k>  «  0.6330  <W/m.<) 

Prandtl  numoer  <Pr>  *  3.834 

Density  <Rho>  •  990.0  (kg/m‘3) 

Kinematic  vis  (Nuv)  »  5.9255E-07  (n'2/s> 

Coef  ther  exD  (Beta)  *  4.1266E-04  ( 1 /V > 

***  Calculations  *** 

Preliminary  calculations: 

Reynolds  number  <Re>  “  7.441E+03 

Friction  ;actor  (stat>  *  3.4Q2E-02 

Heat  flux-  (Goo>  *  4.060E+03  (W/m2> 

Mean  fluid  vel  <Vn>  *  9-23E-9!  (m/s> 

RPM  -  3600.00 

Resul ts : 

Nudb  hunak  hun^  Nuirnn  Mus  Nuwm2 

49.52  55.30  111.80  31.56  44.34  102.34 

Ra  Ro  Ra»Ke«°r  Hf  'lak  1 

7.39E+03  5.16E-01  2^l4E+08  1  .63E-02  44.90 

RaR'e  rg  7v 

5.50E  +  07  7 . 506E-05  0.026 


APPENDIX  C 

THEBBOSTPHON  AHAIYSIS  (COHPOTEB  P BOG RAH  AND  RESULTS) 


The  following  program  was  written  to  give  the  results  o 
the  correlation  for  the  heat  transfer  for  the  rotating 
closed-loop  ther mosy phon  discussed  in  Chapter  III.  It  tn- 
written  in  HP-Basic  for  the  HP-9826  computer.  Sub-programs 
for  the  graphics  and  the  thermophysical  properties  are  no-*- 
shown.  A  sample  run  is  included  immediately  following  th® 
listing. 


Pages  114-117  Program  Listing 
Page  118  Sample  Calculation 


v,  n  z  ri  C_  1 


re rj ii rar*.  -i^- 


COM  I'Csn'  X  r.  i  n  ,  Y  m  >  n  .  3  -  *  .  $  f  •• 
PR  INTER  IS  1 


PRINT  USING  ’ 

■2X. 

.""Default  values:  "" 

** 

D-A.7S3E-3 

» 

Tune  diameter  <m) 

L* .3230 

1 

Tube  length  <m> 

R-.A021 

» 

Radius  of  rotor  <ra> 

Kff-0 

» 

Friction  r actor  »or 

penas 

PRINT~USING  * 

*AX, 

.""Tube  diameter 

""  .  Z  .  ADE  . 

"  (  m  > 

PRINT  USING  1 

’AX 

,""TuDe  iengtn 

"" . DD . D . "" 

( m ) " 

PRINT  USING  1 

•AX 

.""Ckt  Lenath  < l_ t >  * 

•*"  .  DD .  D  . "" 

<  m ) " 

PRINT  USING  ' 
BEEP 

•AX 

.""Radius  of  rotor  - 

""  .  Z . 3D  . "" 

<  n ) " 

INPUT  "OK  TO 

ACCEPT  DEFAULT  VALUES  < 

1  * Y . 0*N  >  P" 

.Id 

BEEP 

INPUT  "LIKE  A  HARD  COPY  < ! *Y .0-N>P" . Ihc 
IF  Ihc* 1  THEN  PRINTER  IS  701 
IF  Id-1  THEN  1280 
BEEP 

INPUT  "ENTER  TUBE  DIAMETER  <m>".D 
BEEP 

INPUT  "ENTER  TUBE  LENGTH  <m>".L 
BEEP 

INPUT  "ENTER  RADIUS  OF  ROTOR  <m)".R 

PRINT  USING  “1  OX . ""*»*  Geometric  Variables  *»* 


1  AX  .  ’ 

"Tuoe 

diameter  (D>  * 

"".Z. 

3DE 

1  AX  .  ‘ 

"Tube 

iengtn  <L> 

""  .DD 

.DD 

1  AX  .  ’ 

"Ckt 

length  <Lt> 

"".DD 

.DD 

1  AX  .  ’ 

"L/D 

a 

""  .50 

.D” 

’  AX  . " 

"Roto 

r  radius  <R>  * 

•  »  *7 

3D. 

PRINT  USING  "  1  AX  ,  ""L/C  -  "".5D.D' 

PRINT  USING  " '  AX . ""Rotor  radius  <R>  *  "”.Z. 3D. 

PRINT 

PRINT  USING  " 1 0X , ""Computed  values  are:  . 

PRINTER  IS  1 
BEEP 

PRINT  USING  "2X, ""SELECT  OPTION: . 

PRINT  USING  "AX, "“I  Single  point . * 

PRINT  USING  "AX. ""2  Nu  versus  RPM . 

INPUT  Op 
IF  Qp>1  THEN 
BEEP 
J;-0 

PRINTER  IS  1 

INPUT  "LIKE  TO  PLOT  U - Y . 0=N > ?" .Okp lot 
BEEP 

INPUT  "Nudb* 1 . Nunak 1*2 . Nuwm*3 . Nus* A . Nul*5" . Ian 

IF  Okplot-i  THEN 

CALL  Piot(Jj.X.Y.Dt.TcJ 

INPUT  "LIKE  A  LABELP(1-Y,0-N)".I1 

IF  11-1  THEN  CALL  Label 

END  IF 

END  IF 

Nstep*1 1 

IF  Op-2  THEN 

BEEP 

INPUT  "ENTER  RPM  RANGE  <  MIN  ,  MAX  > "  .  Romi  .  Rpmn 


*  -  '  ■  ■-  -A  ■ 


kpm-Apmi 


I'^r'jT  ‘  t_  N  T  E  K  nHr1  *  Kt  i  i  npM  4  Pr*'"' 

END  IF 
BEEP 

IF  I he* I  THEN  PRINTER  IS  701 
IF  Ihc-0  THEN  PRINTER  IS  1 
PRINT 

3RINT  USING  “1  OX.  ""*■*»  Operating  Variables  *•» 

PRINT 

PRINT  USING  "  1  OX  , ‘‘“Input  variables  are:  . 

BEEP 

INPUT  “ENTER  HEAT  LOAD  'W)“.Q 

PRINT  USING  "14X.""Heat  ioad  »  "“.4D.DD.""  <H) 

BEEP 

INPUT  “ENTER  COLD  SIDE  TEMP  <C)“.Tc 

PRINT  USING  "1  AX ."“Cold  side  temp  *  ’”* .  4D  .  DD  .  “"  <C>“ 

Dt-5 

Th*Tc+D  t 
Ta*^Tc+Th)»  .5 
Dtt-Dt/L 
PRINT 

Cp*FNCpti<  Ta> 

Mu*FNMuu<  Ta) 

K*FNKui<  Ta) 

Pr*FNPrw<  Ta) 

Rho*FNRhoui<  Ta) 

Nuv*Mu/Rho 
Beta*FNBeta<  Ta) 

Qmega*Rpm*2»PI/S0 
Tn»Tc+Dt 
Ta*Tc+Dt/2 
PRINTER  IS  1 
PRINT  "TA  -  “.Ta 
Rhoh»FNRhow( Th) 

PRINT  "RHOH-" .  Rnoh 
Rhoc*FNRhou<  Tc ) 

PRINT  “RHOC*" .Rhoc 
Rhom*FNRhoiii<  Ta) 

Drho-Rhoc-Rnoh 
IF  Drho<0  THEN  GOTO  3140 
PRINT  "dro*“,Drno 
Mu*FNMuu<  Ta> 

Cp*FNCpui(  Ta) 

Beta-FNBeta( Ta) 

Pr*FNPrw<  Ta) 

Nuv*Mu/Rhom 

PRINT  “MU. CP  -“.tiu.Cp 
Mf-Q/(Cp*Dt) 

Re-4/PI«Mf/(Mu*D) 

Gr*R»(0n>ega*2)*Beta,*<  (D/2  > "  4  >»Dt  t/Nuv"2 

Ra*Gr*Pr 

prod-Ra»Re*Pr 

Nuum* . 262»Prod"  .  173*48/1 1 

Dtwp*Q/<K»Nuwm*pI»L ) 

PRINT  “DTHP  -“ , Dtwp 
DtwDtwp 

Grr*R*0mega  2*Beta*D  3*Dtu/Nuv '2 

Nu*Pr/Gr r*L t/R*,D/L*Rhom/Rhoh»<Kf  f  +  <256/Re*L  t/D )  >*Re' 3 
Dtuc*G)/<PI*L*K«Nu) 


“ :  0 
:  T  c 


116 


I-  ABS<Druc-Dtu:)>.'!p  THEN 
IF  Rpn> 3uG0  'HEN 


OU I U  I Juw 


2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440! 

2450! 

2460 

2470! 

2480! 

2480 

2500! 

2510! 

2520 

a  *• 

2530  * 
2540! 
2550 
2560 
2570 
2580 
2590 
2600 
2610! 
2620! 
2630 
2640 
2650! 
2660! 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
Ta 
2760 
2770 


END  IF 

IF  RpmMOOO  THEN 
Dt-Dt+1 
GOTO  1900 
END  IF 

IF  Rpro>500  THEN 
Dt-Dt+1 .5 
GOTO  1900 
END  IF 

IF  Rpm>-300  THEN 

Dt-Dt+4 

GOTO  1900 

END  IF 

END  IF 

Vm-4«Mf /<Rhom-pI-D '2* 

Ro- Vm/ < Omeqa-D > 

Jay-Omega-D '2/Nuv  !From  STEPH-Rotat  lonai  Re 
Dtt-<Th-Tc>/L 

Gr-R»<Gmega'2)»Beta-<  <D/2>  4)-Dtt/Nuv  "2 
Ra-Gr*Pr 

Laminar  Value 
Nul an-48/1 1 

Turbulent  Value  from  D i t tus-Boe 1  ter 
Nudb- . 023-Re  .  8*Pr' .4 


Stephenson  Correlation 

Nus*<Pr  . 4/ . 367 > * . 007 1  -Re  .38*jay,023  !Pi r . Turb-EQN2 .23-corrected  for  mat 

Nakavama  Correlation  tor  laminar  conditions 
J2-J-2 

Sq-2/1 1+SQR< i+77/4-C 1 /Pr  2) >  > 

Cf-1 - .486-< <3-Sq-1  )  .4/<Sq-<Sq-Pr*SQR<5>-2>>>-J2/(Ra»Re>  .6 
Nunak 1 *48/ 1 1 * . 1 91 /Sq*<  3"Sq- 1 )  ' . 2-<Ra*Re  >  .2-1/<:+(1/1 Q-$q*Pr  >  > 

Mus-Mu  !Temo  Value  Muss'Mu 
Prod-Ra-Re-Pr 

Woods-Horris  Correlation  for  Laminar 
Nuuim-. 262-Prod' .1  73-48/1  1  !  Eqn2.22-HM3 
J1-Jay/8 

Hoods-Morris  Correlation  for  Radial  Pipes 

Nuum2- . 01 5-Re' . 78-U1 ' . 25  !from  wm2 

IF  Ihc-0  THEN 

PRINTER  IS  1 

ELSE 

PRINTER  IS  701 
END  IF 

IF  Jj-0  THEN 

IF  Jj-0  THEN  PRINT  USING  "1  OX . ""Resul ts :  . 

PRINT  USING  "1  OX , ”"F iuid  properties  evaluated  at  ””.DDD.DD,"”  <C>  are: . : 


'.DDD.DD.""  <C>  are: 


°RINT  USING  "  1  4X  . ‘‘"Spec i f  ic  heat  <Cp> 
PRINT  USING  “1 4X  . ""V iscos i  ty  <Mu> 


"" , 4D . DE . ( J/k  q . K ) . : Cp 

""  .  Z  .  4DE  .  ""  (N.s/m  2> . :Mu 


'  I£tr- 

-PINT  "3? IMG  " 

i  sX  ,  ""T'lerna :  ;cnn  <  t  > 

-  "".Z.JLV"  <«/n.p  1 . . 

27gri 

PPJNr  USING  '• 

1  4X  , ""Pranat  i  numoer  <Dr> 

*  ""  .  Z  .  3DE"  : p r 

:r-.oo 

■  2 S 1 1  i G 

1 4X . ""Dens i ty  (Rho> 

-  .  4D  .  2  .  9k  a/  n  2  / .  . 

2  _  *  v 

1  4X  .'"'Kinemat  ..p  (Nu..  1 

J.OC.U 

-•'*.14  1  Lp.lNb 

MX,""Uoef  therm  exp(De:a> 

*  ,  p  .  upp  .  •••  <  i  /«.  /  :  pe 

2830 

PRINT  USING  " 

1 4X , ""Reynolds  number  <Re> 

-  .  Z  .  3DE"  :  Re 

2840 

PRINT  USING  " 

14X."”*1ean  fluid  vei  <Vm) 

-  ""  .  Z  .2DE  .  ""  (m/s) . •2-M 

2850 

PRINT  USING  " 

i4X.""Mass  flow  rate 

-  .  Z  .2DE  (ka/s> . :M 

2860 

PRINT  USING  " 

14X.  Temperature  <  hotsioe)  -  ""  .DDD.DD  .  ;  :  h 

2870 

PRINT  USING  - 

14X,""RP*i 

-  ,  DDDDD  .  .  "  :  Rom 

2880 

END  IF 

2880 

PRINT 

2900 

IF  2kpiot-0  OR  (Okpiot*l  AND  Ns*.ep<11) 

T'-EN 

2810 

PRINT  USING  " 

10X,""  Nudb  Nunak 1 

Nuum  Nus 

2920 

PRINT  USING  " 

*  3X , 6  <  3D . DD , 2X)" : Nudb , Nunak 1 . Nuwm . Nus 

2930 

PRINT  USING  " 

1 4X , ""  Ra  Ro 

Ra-Re-Pr  Mr 

2840 

PRINT  USING  " 

i4X,n(Z.2D£.2X) .2D.DD":Ra. 

Ro  .Prod . Mf 

2950 

END  IF 

2960 

IF  Okpiot*1  7 

HEN 

2970 

IF  Ian-1  THEN 

Y-Nudb 

2980 

IF  Ian-2  THEN 

Y- Nunak 1 

2990 

IF  Ian-3  THEN 

Y-Nuum 

3000 

IF  Ian=4  THEN 

Y-Nus 

3010 

IF  Ian-5  THEN 

Y  =  <  48/1 1  ) 

3020 

X-Rpm 

3030 

Jj-Jj+1 

3040 

CALL  PIot<Jj .X, Y.Dt.Tc) 

3050 

PRINTER  IS  1 

3060 

END  IF 

3070 

IF  0o-2  THEM 

3080 

Dt-5 

3090 

Rpm-Rom-I 0.1 

2100 

IF  Rpm>Rpmh  7 

HEN  3130 

31 1  0 

GOTO  1890 

3 '20 

END  IF 

V  30 

BEEP 

3140 

INPUT  "ANOTHER  RUN  < 1  * 7 . 0*N >7" . Ir 

3150 

PRINT  "PU" 

3 1 50 

IF  Ir-1  THEN 

1360 

3170 

INPUT  "WANT  TO  LABELS  1  -Y  .  0-N)'M  I 

3180 

IF  11-1  THEN 

CALL  Label 

3190 

END 

a  ho 

ra 


:  •  :  "  j  r  ir 


TT7TTT?  777 


iV 


i 


i 

,v 

V 

N 

> 


i 

s' 


[  * 
r; 

fc- 


.‘c  f  *C  /  a. .  .  v^.  *  c* 

Tone  diameter  (D>  ■ 
f  ij.;,o  length  ( t_  >  * 

Ck  t  length  (Lt) 
i_/D 

Rotor  radius  <R> 

Computed  values  are: 

•»■*  Operating  'lari ad.es 

Input  variables  are: 

.-eat  .  o  ad 
Cola  side  temp 

Results: 

P’luid  properties  evaluated 
Specific  neat  (Cp> 
Viscosity  (Mu' 

Thermal  cond  <k  > 
Prandtl  number  (?r> 
Density  (Rno> 

Kinematic  vis  <huv> 
Coef  therm  exp(Beta) 
Reynolds  number  (Re) 
Mean  fluid  vel  (Vm) 
Mass  flow  rate 
Temperature  <hotside> 
RPM 

Nudb  hu.nai  1  huwn 

2.42  51.30  34.41 


4 .  /b3E-0?  (m> 
.92  (n> 

2.45  <ra> 
172.3 
0.402  (m > 


-  SO.flC  («> 
45.00  (C> 


at  £2.75  < C >  are: 

-  4310.1 E+00<  J/k° .K> 

=  4.41 72E-04  <  M , s/m  2> 

-  0 .6404  <W/m.k> 

-  2.902E+00 

-  939.1  9k  o/r.  3> 

=  4.4991E-07  ( M  ' 2 / S > 

*  5 . 3940E-04  (1/K) 

-  1  . 978E+02 

-  1  .87E-02  (m/s) 

*  3.27E-04  ( k a/s ) 

-  80.50(0 

-  3600 . 

Nus 
1  .65 


Ra  Ro  Ra*Re*r'r  Mf 

6 . 1 3E+05  1.04E-02  3.52E+08  3.2 ?E-0« 


*1 


t 


w 

3 


a 


■i 

n 
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A P PEP PIT  t> 

CLOSED  TUO-PHASE  T H2 S30SYPH  OH  ANALISIS  (COMPOTES  PROGRAM  AND 

RESULTS) 

The  following  program  was  written  to  aive  the  results  of 
the  correlations  for  th«  sonic,  boiling,  and  condensincr 
limits  as  discussed  in  Chapter  IV.  It  was  written  in 
HP-Basic  for  the  HP-9826  computer.  Sub-programs  for  the 
graphics  and  the  t hermophys ical  properties  are  not  shown.  A 
sample  run  is  included  immediately  following  the  listing. 
-Pages  120-124  Program  Listing 
-Pages  125-127  sample  Calculation 


1200  BEE3 

'2’0  I'lP'j-  ‘".IKE  A  -ARD  COP7  <  '  =v . 0»N> 3" . Ihc 
1220  I3  I  nc-'  ~HEN 
■230  PRINTER  IS  TO! 

1 2 AQ  ELSE 
r 250  PRINTER  IS  ; 

1260  E! ID  IE 


1270 

IF  la- 

1  THEM 

1360 

i  25 C 

BEE3 

1220 

input 

■'ENTER 

TiJEE 

DIAMETER  < . D 

1300 

9Er- 

1310 

INPUT 

"ENTER 

TUBE 

L£MGTH<£  >ap >  . 

1320 

BEE3 

’3  30 

T>|P»JT 

"ENTER 

TUBE 

_ENGTH<Cona>  <n)”. 

£ 

i  340 

SEE3 

:  :so 

INP'J  T 

"ENTER 

RADII 

IS  3F  ROTGR  R 

'360 

PRINT 

USING  " 

•  1  OX  .  • 

'"»*■*  Geometric  Var 

iacl.es  ••• 

II  M  M 

:  370 

PRINT 

USING  " 

1  1 4X  .  " 

'"Tiioe  diameter  <0> 

-  "" . Z . 3DE 

.""  (  m )  ’" 

;  D 

1330 

PRINT 

USING  " 

’14X." 

'"IiiDe  ieng  t  Ti  <  Evao  > 

<Le>  -  "" 

. Z . 3D . "" 

<  m  > 

1  390 

3RIN? 

USING  " 

'  1  4X  ,  ""TijLDe  lenatn(Cona) 

<Lc)  -  "" 

.  Z . 3D . "" 

( n  > 

1400 

PRINT 

USING  ’ 

*  1  4X 

'"Rotor  raGms  <R> 

_  it  •• 

. Z.3D."" 

( m ) 

141Q 

PRINT 

1  420 

PRINTER  IS  1 

1430 

BEE3 

i  440 

3RINT 

USING  " 

’2X . "" 

'SELECT  OPTION: . 

1450 

PRINT 

USING  ' 

•4X 

Q  versus  /.Fill"' 

T  460 

PRINT 

USING  " 

'4X.""2  R??  . 

i  470 

INPUT 

0p 

1  480 

IF  0o> 1  THEN 

1400  BEE3 
1 500  PRINTER  IS  t 

1510  PRINT  USING  "4X. 10 T  IMPLEMENTED  YEr~ 

'529  GOTO  14Zi!) 

1530  ELSE 
1540  BEE3 
1550  G-9 

1 550  J  '*0 
'570 

1530  INP'J T  "LIKE  '0  3L0T  (l.y.O-i)T".OkDl3t 


O  19 


. »0.:< .  1 .  rype.Cx> 


20 

3PPP 

:  0 

INPUT 

,40 

END  IF 

*50 

END  IF 

eQRSH~  3P':>H  .  3  s 


N-31-O0*  TOO 
5EE- 


'630 

TNPUT  "ENTER  OP  TON  <  !  =S 

,2»E.3*3.p»C 

1690 

Ic  In-2  THE*: 

t  700 

BEEP 

’7*0 

INPUT  “ ’ -HWllIS . 2= J ASTER’ 

"  .  le 

<  /  jl 

END  IF 

:  ’TO 

"  THEN  PRINTER  IS 

70 1 

1  7  40 

it  iHc=n  turn  PRINTER  IS 

* 

’’50 

BEE3 

:  75 0 

:.NP<!7  "cNTER  'sat  <C>"  ■  ~ 

_  - 

t 

1  7’0 

PRINT  USING  "14X.""Tsat< 

c  > 

'  730 

BEE3 

17  90 

INPUT  "ENTER  RPM" . Rpm 

’  300 

PRINT  USING  ”*  4X , ""RPM 

1310 

3EEF 

TEMPERATURE*1 

’  320 

INPUT  "ENTER  TEND-  -4ALu 

1  330 

PRINT  USING  "*  4X . ”"T  mal 

l 

:  340 

At*°I-0  2/4 

1350 

Tsatk  =  1  sat *273 .  '  5 

'  360 

2p*FNCduI 'sat ) 

1370 

Mai  ** NMaw<  <  T-sat  +  Tw) /2  * 

'330 

11uv-“NMuv<  Tsat) 

'990 

K-rNK«i  (  T sat  > 

1900 

Pr**- NPrw< Tsat ) 

19'0 

Rhoi=FNRhow<Tsat) 

1920 

rlf  g»*NHf  g<  Tsat  >  ’ 

1930 

3eta=FNBeta<  Tsat  > 

•  94Q 

S igma=*NS :gma( Tsat  > 

1 950 

Rhoy -  * /FNVvs t  <  Tsat ) 

I960 

Mav “May/Rhow 

1970 

Rsn* 1 .327 

1980 

Cw'Co/Rsn 

'990 

Rgc*4Bi.52  1  <  J  /  k  g  .  K  > 

2000 

0nega*Rpm*2*PI/6G 

2010 

Acce i -Omega  2*R 

2020 

IF  J  j  *  0  THEN 

2930 

PRIM' 

2040 

PRINT  USING  "10X.  Fluid  orooerties 

’ .  DDD  .  D‘ 


.  -<D" :  Ren 


’ .  DDD  .  D’ 


\DDD.DD. 


(C>  are:' 


"  •  3  <3.  t 

E050  =>ri: 

2060  PRI 
2070  PRI 
2030  PRI 
2030  °RI 
2iQ0  PRI 
2’ t0  PRI 
2120  PRI 
2 ' 30  PRI 
2160  PRI 
2 '50  PRI 
2160  END 
2' 70* 

’•'■or  ?<"- 


PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
END  IF 


’ 1  4X  ,  ""Soec if ic  neat  <Ce> 

‘14X. ""Ratio  of  Specific  Heats* 
*1  4X  .••■'Viscosity  <Hu> 

’1 4X . ""Thermal  cone  <k> 

M  4X .  ""Prandti  nunoer  <Pr> 

*1 4X , ""Entnaloy  <rifg) 

’  i  4X  .  ""Dens  i  ty  <Rhoi> 

M4X. ""Density  <Rhov) 
*!4X.“"Coet  trier  exo  'Beta) 
*14X. ""Omega 


\4D  D."M  (J/lca.K) . :Cd 

’ .  D . DDD" : Ksn 

‘ .  Z  .  4DE  .  (M.s/m  2> . ’.Mai 

' .  Z  .  4D  .  ""  (U/m.K) . :K 

’ .  Z  .  3D” : Pr 
’ .  Z  .  4DE” ; Hf a 

' .  4D . D . ""  <’<g/n  3> . : Rhoi 

*.Z.40.‘"'  <  It  a/m  3> . :Rhov 

*  .Z . 4DE  .*•”  ( f  /K) . -.Beta 

1 . 4D  .  D . ""  <i/s> . : Omega 

‘ .  Z  .  4DE  . <N/m> . :Stama 


*  ‘  ‘-/iL*  ~  ‘/-'iL*.  -  *■  ».*■  k*>I  ilV  - 


22  o  o 

2  n 

2220  • 


r>  Kov-r'  "o  . )  )  )  *■  <  t  Rno  ;  •  ^  f’ ho*^“^^o  •.  ^  )  *  (A  t ) 


2230! 

SONIC  CALCULATION 

2240 

IF  Io= "  !HE: 

2250 

IF  2--0  THEN 

2260 

PRINT 

p  ^  7  0 

°RTNT  jSING 

•••OX  .“"Srtiic 

2230 

PRINT  jsiNG 

"i ox. -••/:  r i i 

2290 

END  " 

2300 

K  '*  /  ( 1 1  *  <  *’h<W 

23  iq 

MT -Om/Ht 

2320 

X-p.;- 

2230 

Y  -On 

2340 

IF  :3n<»9  THEN  GOTO  2400 

2350 

IF  Okoiot*1 

THE!-; 

2360 

CAL_  ?  I o  t  <  J ; 

•  Jc.-a.X.Y.  r 

2370 

ELSE 

2330 

PR  IN T  USING 

" f  2.X  .  3D  .  3X  .  Z 

2390 

END  IF 

2  400 

IF  °cf >= ' 00 

THEN  GOTO  37 

2410 

IF  Okoiot-1 

THEN 

2420 

*cf * i 0  . 

02 

2430 

ELSE 

2440 

P c t  * p c *  *  ? 

2450 

J ; *  j ;  + 1 

2460 

END  :=■ 

2470 

GOTO  2200 

2430 

END  IF 

Total  He= 


2^*30  ? 

2500 !  ENTRAINMENT  CALCULATIONS 

2510  IF  Io-2  THEN 

2520! 


2530 

2540 

2550 

2560 

2570 

2580! 

2590 

2600 

2610 

2620 

2630 


2640! 
2650! 
2660 
2570 
2680 
2C30 ! 
2700 
2710 
2720 
2730 
2740! 
2750! 
2750 
•>770 
2730 
27°0 


IF  J;-0  THEN 
PRIMT 

PR  INI  USING  "1  OX  .  ""Entraining  Limit  Claculations  :* 

PRINT  USING  Fill  Total  Heat . 

END  IF 

HALLIS  CORRELATION  >-.25 

IF  le-i  THEN 

IF  AvO  THEN  GOTO  2990 

Jv-  . 25/Rhov  .5*< Omega" 2*R*D*(Rho:-Rhov))  ' .5 
Gm-Rhov*  Jv»At»Hfg 
END  IF 


JASTER  CORRELATION 
IF  Ie-2  THEN 
IF  Pot-0  THEN  Pc  f -  I 
T1 -< 1 /( I -<Pcf /100> )  )-i 
T1-< ! / < Pcf / 1 00 > » - ! 

T2-T1 »<  Rho  1/Rhov ) " ( 2/3 ) 

IF  T1<0  THEN  2990 

Xx2*<Mui/Muv>  ’  .25*1'  •  .75*<Rhol/Rhov>  <  < 2/3 ) - 1  .75-1  > 
Xx-i /( 1 +72) 

Xx»  . 5 

<x2  -  <  Mu.1  /Muv  j  ' .  25*Rhov/Rho  1 

IF  G-0  THEN  G-i 

Re  1  -  <  !  -Xx  >  ■*<  G*D/Mui  > 

Rev-Xx-b-D/Huv 
Fy  .  0  79/Rev  .25 


7 W/ 


2':'  : 0  IF  -\e-M_5U  ’HE.'I  jO3  ■'  ..; 

2820  Nun*  “  A-»FG^  {  r  >■*  t  Rh*'-  i  »Pfi;v  > 

2930  >rs«<  ?  «-<x2*  2C-  >  -s.S-hnX.  •  : 

25  4  U  .J  C  *  1  N  ill !>  ‘  J »?  71  *  I  .  .*  I  J 

21-50  ir  ABu  i  v  ’.i  -ijO  ^  -  be  ’ '  •  j  i  nE.'l 

2360  G=  < b  mjc > /_ 

2370  G0T0  2760 
2830  END  IF 
2890  Qm-Xx»G-<U-Hfg 
2900  END  lc 
29 :0  <-Pcf 
2920  Y-On 

2930  IF  Qpk-0  THEN  GOTO  2990 
2940  IF  Okoiot*"  THEN 

2950  CALL  3lot  <  j  ;  ...ic  .  .  ■< .  Y  .  T/pe  ,Cx  > 

,?<3|j0  r!  ''c 

29’9  PRINT  USING  ”12X  .  3D.3X , Z .4DE ,4X . Z .3DE" :?cf .Qm .0 
2980  END  IF 

2990  IF  Pcf  >* 1 00  THEM  GOTO  3720 

3000  IF  0k oior'*'  THEN 

30 10  Prf  =  s.;+*!  O  ' .02 
3020  ELSE 

J  0  ^  0  ?r*' ■*<<•+  +  !  0 

3040  END  1^ 

3050  J  r -  : 

■060  IF  Ie*2  ’HEM  GOTO  2680 
3070  IF  le-1  THEN  GOTO  2200 
SOSO  END  IF 
3090! 

3  ’  00 !  SOILING  CALCULATIONS 
3110  IF  Io*  3  THEfi 
3 '20  IF  J;-0  THEM 
3 1 30  PRIME 

3 !  41)  PRINT  USING  "1  OX .  "”3o  1 1  ing  Limit  Caictiiac  ions : 
3150  PRINT  USING  "■OX."-/.  Fill  Toral  Hear  QDP 

3160  END  IF 

3170  IF  °c: > ‘ 00  ’HEM  GOTO  3720 
3180  Ab*D*'_e*ACS  (  1  -Pcf /50 ) 


imiii 


3’90  Qm-Ab- .  I  3»Rhov/  .  5*Hf  g-XQmega  2-*R*<  Rho  i -Rhov  )*S  j.gma>  '  .  25 
3200  X-Ocf 
3210  Y-Qm 


3220  IF  Qm<  *0  THEM  3230 

3230  IF  0kDiot*1  THEN 

3240  CALL  3  lot ( J j . Jc . Jd . X . Y . Type ,Cx > 

3250  ELSE 

3260  PRINT  USING  " 1 2X . 30 . 3X . Z . 4DE ,3X . Z . 4DE" : Pcf  .  Om .Om/Ab 
3270  END  IF 

3290  IF  Pcf >  1 00  THEM  GOTO  3730 
3290 

3300  IF  Okoiot-l  THEM 
3310  Pcf “°cf *i 0  .02 

3320  ELSE 
3330  Pcf -Pcf *1 
3340  END  IF 
3350  GOTO  2200 
3360  END  IF 
3370  * 

3380!  CONDENSATION  CALCULATIONS 
3390  IF  Io-4  then 


124 


34 ’ 0  PR  IN 

3420  DP;MT  (.'SING  OX  .  ‘"‘Concensi  na  Limr  Calculations 

3430  PRINT  USING  "’OX.'""/,  -i.il  local  Hear . 

3440  END  IF 
3450  Dt=!sat-Tw 

3460  HF  gp-Hf  1  *■ .  5S*<  Cp-Dt /Hf  g  >  > 

3470  Tag*  <  Rho  i * <  Rho  1  -Rhov  >  •Qmeoa  E^R-Hf  ap*i'  3/ ( D*Mui  * 
3400  DEG 

3490  ?Hi*'30-AC3<  !  -<°cf  /50>  > 

3500  IF  Phi>*0  AND  Pni<-H0  THEN  F*CNF  i  <Ph  i  > 

3550  IF  ?Si>l 10  AMD  P  h  i  <  = ! 5  0  THEN  F*FNF2<?hi> 

3520  IF  Phi>150  AND  Phi<*’,30  THEM  F»-NF3<=hi) 

3530  IF  °hi>5S0  'HEN  PRINT  "ERROR .PHI> ! 30 ! ! ! " 

3540  RAD 

.3550  Qn*Av»Dt*F*Toa 
3560  X-Pcf 
3570  v=i3ni 

3580  IF  Qm<*0  THEN  GOTO  3700 

3530  IF  Qkpiot*1  THEN 

3600  (3ALL  P Lot < U j .Jc  .  Jd .X  .  Y . Type  .Cx > 

3610  .  ELSE 

3620  0RI?IT  USING  " • 2X . 2D . 3X . 3D .  D" : Per  . On 

3530  END  IF 

3640  IF  Oko lot* 1  rHEM 

3650  Dc^*Dc*'*-0' .02 

3660  ELSE 

3670  Pc  * *pcf  + 1 

3680  END  IF 

3630  J  ’ *.J  r  + 1 

3700  IF  t > * ’ 0 0  THEN  GOTO  3730 
3710  GCTO  2200 
3720  END  IF 
3730  BEE? 

3740  INPUT  "ANOTHER  RUN  < ! - Y . 0-N > ?". Ir 
3750  IF  Ir-i  THEM  1200 
3760  BEE? 

3770  INPUT  "'4ANT  TO  :_ABEL?<  1  *  Y  .  0-N)  ".I  I 
3780  IF  Ii=1  THEN  CALL  Laoei 
3790  END 
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